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SYITOPSIS 

The present thesis is divided in four chapters. 
Chapters II a.nd III describe the major contribution of the 
present work. 

Chapter I extensively reviews the work done during 
the last two decades in approximation and application of 
Case’s theory in standai’d neutron transport problems. In 
the li^t of the review, it tries to establish the objective 
of the present work. 

Chapter II develops the formulation of Transport 
theoretic Pjj approximation (TP^), The methodology is based 
on the framework of conventional Pj^j-approximation in conjunc 
tion with Case’s exact treatment of the one-speed neutron 
transport equation. Firstly, it demonstrates that the use 
of all the roots of conventional Pjj approximation is not 
mathematically justified and the exact asymptotic decay 
constant must be used in place of its counterpart 
in describing the asymptotic flux. Formal solution of the 
transport equation is then constructed with the aS 3 niiptotic 
part being of the same nature as exact solution and the 
transient part terminated as in conventional Pjj approxi- 
mation. The second part of the scheme generalizes the 
liarshak’s free-surface boundary condition. This is achieved 
by introduction of a parameter in the fi3?st equation of 



the set of Ifershak boundary conditions, a^j- is evaluated 
by forcing the asymptotic flux for the source-free Milne 
problem to be identically same as the excict transport 
theory asymptotic flux. The Milne problem is chosen as a 
reference problem 3 ,s it leads to the free surface condi- 
tions. The same can be used for any other dne-speed 
transport problem involving a free surface. Analytical 
expressions for the parameter for TP^^ and TP^ approxima- 
tions are presented, hhile analytical expressions for 
can be similarly constructed for higher order approximations, 
a simple formula is prescribed for evaluating (F ^ 3) to 
avoid prohibitively-increasing algebraic complexity. The 
formalism is then applied to different standard transport 
problems - the emergent distribution and leakage for the 
Milne problem, the scalar flux in Milne problem when c = 1, 
the leakage in constant source Milne problem and the criti- 
cal problem. The calculation shows that both the emergent 
distribution and the leakage for the source-free Milne 
problem improve significantly in the TPjj approximation. As 
a check on I the justification of the present scheme for the 
special case of c =1, the scalar flux obtained by this 
scheme for Milne problem is compared with the exact and 
Pjj values. TP^ flux is superior to P^^ flux for large 
values of X, while TP^ improves P^ calculation for all x 
except for very small values of x. Constant source Milne 
problem and th® erf. 'tica-lity problem are next considered 
where, unlike the source-free Milne problem, a^j. is explicitly 



used in th.e formulation. In the constant source Milne 
problem it is found that method worsens the leakage 
value for low values of c. The reason for the failure of 
TP|,j method in this C3,se may be attributed to poor tra.nsient 
flux representation and also to the fact that extrapolation 
distance for constant source problem for small values of c 
is widely different from that of source-free Milne problem. 
As is apparent from calculation, TPj^ significantly improves 
the critical slab thickness. TP^ gives the end-point 
result and TP^ improves the result further. TP^j- theory is 
then applied to modify the interface conditions and a typi- 
cal fuel-moderator cell problem is then analyzed. It is 
found that TPjj-approximation significantly improves the 
disadvantage factor compared to its Pj^ counterpart, Finally 
it is concluded that the present scheme usually leads to 
improved results except in cases where the transient compo- 
nent is overwhelmingl3;' predominant. Also, unlike asymptotic 
diffusion theory, any order of approximation can be syste- 
matically constructed in this method. 

Chapter III poses the finite slab problem in a multi- 
plying medium as an inverse eigenvalue problem and treats 
it functional analytically to establish the existence of 
the critical thickness. Numerical estimate of critical 
thickness is also obtained by functional analyiiic methods. 
The set of two coupled Fredholm equations as obtained in 
Case’s theory for a finite slab problem with an incident 
beam on one of the surfaces is combined together to form a 



single Fredholm integral equation* Critical slab problem 
is defined by putting the non-homogeneous term equal to 
zero. The equation is then recast in a simpler form for 
mathematical analysis. This is analyzed as a general eigen- 
value equation in which an operator IK('b) depends on a 
parameter b. The problem we are interested in is that one 
which specifies the eigenvalue and determines the parameter 
b (i.e, the operator). Hence it is caHled the inverse 
eigenvalue problem. The operator completely con- 

tinuous in except for those b where the kernel has 
singularity. , 33 )" straight-forward application of Gohberg 
theorem in the range of anal 3 rticity of the operator-valued 
function \K(^)s "the existence of critical thickness is 
proved, Analyticity criterion on lK,Cb) shows that the 
critical thickness b^ must be less than the end-point 
result. This fact is not reflected by Pjjj. calculation. The 
kernel K(v , b) of the integral equation ( v, e(0,l) 

are dependent variables of the integral equation) is then 
analyzed for positivity in for all b. Restriction 

on positivity is necessary for elegant application of 
bounds on spectral radius of linear operators. It is 
established that an approximate problem which sli^tly 
departs from the original equation (numerical calculation 
substantiates this) satisfies the positivity criterion. In 
the general problem, ^(b) = 1 defines the critical slab 
problem and our intention is to solve the equation Kt) = 1 
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for critical half -thickness b^. Estimates of spectral 

/ 

radius in its simplest form are used to find the critical 
thickness. Considering the crudity of the bounds of spec- 
tral radius, the result is found to be fairly accurate. 

Ko sophistication in calculational technique has been 
attempted to improve the result fiirther. This is because 
that the present work is not meant for the calculation of 
critical thickness but is aimed at precise foimulation of 
a class of inverse eigenvalue problems. Numerical calcula- 
tion is only done as a check on the hypothesis adapted in 
the proposed analysis. The choice of the critical problem 
is only incidental and the theory can be applied to the 
general class of similar problems. At the end of the 
chapter, a discussion oh the existence of solution for a 
slab albedo problem is given. 

In the last chapter, a separate Transport-theoretic 
Pjpapproximation scheme has been constructed for constant 
source type problems. It has been shown that the escape 
pro1)ability results for a finite slab containing uniform 
source are significantly improved in the present method 
over conventional Pjj-approximation, whereas TPjj-formalism 
based on source-free Milne problem fails absolutely in 
this case for low values of c. A discussion on possible 
further extension of the present work is also given. ' 



CHAPTER I 


nTTRODUCTIOH 

The pherLomeiaon of neutron transport, which, is of 
primary significance in the analysis and design of nuclear 
reactors, is adeq_uately described by the linear Boltzmann 
equation. In the most general form of the Boltzmann equation 
the neutron angular density appears as a function of space, 
angle, energy and time. Unfortunately this equation in 
its entire generality is not solvable analytical].y in most 
cases. Thus, one has to consider the problem in a restricted 
subspace with drastic simplifying assumptions so as to be 
able to treat neutron diffusion theoretically. The present 
thesis restricts the analysis to one speed neutron trans- 
port in slab geometry. In this approximation, the neutron 
speed is assumed to be unchanged in a collision though 
angular deflection occurs. The assumption may seem to be 
severe to describe the physics of neutron transport but 
it is well-known that 'the mathematical analyst is no 
respector of physics if lie sees a chance to obtain an 
exact solution' .[l] Nevertheless, one-speed theory has 
great importance for various reasons. Firstly, analytical 
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Glosed-form solutions can be obtained for different 
standard problems. These solutions provide a thorough 
mathematical understs^nding of the Boltzmann equation s-nd 
is exceedingly useful as a qualitative guide for practical 
reactor calculations. The standard problems also provide 
a sample set for testing computer codes. The more import- 
ant factor is that the solution of one-speed equation 
forms the basis of multigroup approximation, frequently 
used in the design of nuclear reactors. 

The framework of the present thesis is based on 
exact and approximate methods of solving one-speed neutron 
transport equation. For exact treatment, we choose 
Case's singular eigenfimction expansion technique while 
spherical harmonics method is our choice for approximate 
analysis. A great deal of work has been done on Case's 
singular eigenfunction method and the technique has 
become almost standard, in preference to the older 
VJiener-Hopf technique in the exact treatment of neutron 
transport. From the point of view of analysis, the 
Pjj-approximation is possibly the most popular approximate 
method of one-speed transport equation solution. As is 
well-known, the Pjj- approximation tends to exact solution 
as FT -► , and the ob jective behind improvements in 

Pjj-approximation is only to achieve faster convergence i.e, , 
more accurate solution with a lower order of approximation. 
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During last two decades , much effort has been made to 
improve the Pjj-approximation by incorporating transport- 
theoretic effects in it. However, in these works better 
angular flux representation through modification of usual 
boundary conditions in a general way did not receive 
proper attention. 

The present work is broadly divided into two parts. 
In Chapter II, a Transport-theoretic Pjjj.-approximation 
scheme has been formulated for the approximate solution 
of one-speed equation. The asymptotic part of the solution 
has an exact angular flux representation while the tran- 
sient is still in the framework of Pu-approximation. 
Numerical calculations of different standard transport 
problems have been performed. Through a better angular 
flux representation, faster convergence of solution is 
achieved. In Chapter III, the critical and slab albedo 
problem have been analysed by using Case's method and 
function analytical theorems. The critical problem is 
studied as an inverse eigenvalue problem of a Predholm 
integral equation. Numerical calculation of critical 
thickness has also been performed by function analytical 
methods. The results are good considering the simplicity 
of the method which does not involve the actual solution 
of governing Fredholm equations. 
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Lj-uera turo Revie w; 

Tho literature roviou on exact solution of one- 
specid transport equation is ^iven first to enable a better 
imderstanding of approximate metlaods in the light of exact 
analysis. 

Among the early xforks in the exact treatment of 
one-speed transport equation, the first noteworthy technique 
was developed by Wiener and Hopf [2], In this method, a 
functional equation is obtained by suitable Fourier or 
Laplace transform of the original one-speed equation and 
this equation is then solved by the method of complex 
analysis. The method was first successfully applied for 
the solution of Milne problem in the case of radiative heat 
transfer and then for different boundary value problems of 
neutron transport theory. However, with the advent of 
Case’s work [3] on singula-r eigenfunctions, the Wiener- 
Hopf method lost its supremacy. The idea of singular 
eigenfunction expansion technique wa.s first put forward 
by Davison [4] and then elegantly developed with proper 
mathematical sophistication by Case, The basic idea of 
the method is to separate the space and angle variables 
by assuming tho angular flux 'f' (x, iJ.) of the form 
’J’(x, p) = and then to construct a complete 

set of eigenmodes in a suitably chosen function space. It is 
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demonstra-ted that the spectrum of v is a set of two discrete 
eigenvalues + V Q outside the interval (-1; 1) and a continu- 
ous part o (-Ijl). The formal solution of the transport 

eq.uation is then constructed as 

-x/v x/v 

!>) = ^ 0 + • ® ® 


+ f A(v ) e ^ (p) d V 

'-1 

where (p-)i. ‘{“v (P-) ^-re the eigenfunctions corresponding 

to discrete and continuous spectra of v respectively. 
and AC'^) are the expansion coefficients which are to he 
determined from source or boundary conditions. For a full- 
range problem, this condition is of the form 


f (p) 


a 


O-f ^04- 


O-i) + 4>o- (p) + j” 'f'v(p) dv 


The full-range completeness theorem proves that the functions 
*^ 0 - 1 - are complete for a given function f(p) 

defined on the full range of p e (-1, 1) and determines 
the expansion coefficients by orthogonality relations. [5] 

For a more realistic problem consisting of a half-space 
only, the boundary condition takes the geneiral form 

f(p) = ({)q^ (p) + y A(v) 4)y (p) dv 0 1 W <1 
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Th.e expansion function fCp.) defines a particular problera 
as for example, f(p) is 6(!i. - for half-space albedo 
problem and - in the Milne problem. It can be shown 

by half range completeness theorem that the functions 

^ L I. i 3-^® complete for the expansion 
function f(iJ.) defined on the range 0 p ^ 1. The expe,nsion 
coefficients are determined by application of the associated 
orthogonality relations [5]. This basic methodology of 
Case's work is discussed elaborately in the texts of Case 
and Zweifel [5] and Bell and Classtone [6]. Applications 
of Case's method are given in r ef erence- [ 5] and a review 
article by McCormick and Kuscer entitled 'Singular eigen- 
function expansion technique' [7]. Case's formalism for 
isotropic scattering is extended to the case of anisotropic 
scattering by Mika [8], Mendelson and Summerfield [9] 
applied the singular function method to problems involving 
two adjacent half-spaces. They succeeded in establishing 
suitable orthogonality relations similar to those obtained 
for a single medium. The completeness of eigenfunctions 
for multigroup neutron transport was taken up by various 
workers [10-14]- Numerical calculations oh multigroup 
transport problems were done by Pahor and Shultis [10], 
Siewert and Ishiguro [11], Metcalf and Zweifel [12]. The 
computations, in general, are much involved and it is 
doubtful that singular eigenfunction expansion method can 
be elegantly applied for worthwhile numerical calculations , 
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The discussion presented in the earlier section is 
"based on either infinite medium or half- space problems. 

As the present work is concerned with the finite slab, we 
emphasize this pa^rticular problem. In his formulation, 
IhLtsis [ 15 ] considered a finite slab of half-thickness b 
in a multiplying medium (i.e., c, number of secondaries 
per primary, greater than 1) a-nd derived a singular integral 
equation in terms of the expansion coefficients using the 
symmetry and boimdary conditions of the problem. The 
singular integral equation was then transformed to a set 
of two coupled Fredholm integral equa.tions by standard 
methods [16], The solution of the problem is given by the 
following equation with a^^ = a^_ = 1, say, 


^(x, p.) = (x, p) + (x, p) 


+ . A(v ) iPy (x, p) 

'^-1 


A('^) obeys the following Fredholm integral equations 

b/p b/v —2b/ .y 

e A (p) = -e ° ^ ^ (-v_)e °] 


2 ; (vq^- p^).(l-c) X (-p) g(c,p) - ( v^^-p^)(l-c) 

X gCc.ti) z (-^)/ a, 

(v + p) 


( 1 . 1 ) 
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V e 
0 


[2 


(Vq) - O 


-2b/^ 


X (. 


-o)] 



a(v) 


e- 


V X (-v) dv 


( 1 . 2 ) 


where, X(v) is the well-kiaown X-function of Case's theory and 
g(c, li) is given Toy 

^ ^ 1 

( 1 - cp- tanh“^ pi ) ^ + -2.^-- 


Equation (1.2) is the exact criticality relation - a relation 
Toetween the material concentration (c or v^) and the 
thickness (2b) of the slab. Thus A(v) has to obey the 
Fredholm equation (1,1) and at the same time satisfy the 
criticality condition (1,2), Mitsis constructed Neumann 
series solution for equations (l.l) and (1,2)* However, 
because of the complexity of the nature of equations, the 
iterations could not bo extended beyond first order approxi- 
mation (which itself is rather complicated). The equation 
(1,1) is written in the form 


A (v) = A, ( 


1 

v) + X 
• Jo 


K (v, a) A (a) da 


the Neumann series soltition of which is given by 

.1 

'o 


A (v) = A^ (v) +x( K( V, a) A^ (a) da 

-^o 
1 ^1 

\ ] K (v, a) K (a, p) A, (p) da dp 

^o-'o 


+ x*- 

+ # • * 


(1,3) 
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S’or a zero-th order or diffusion-like approximation where 
there is no tra,nsient component, Mitsis assumed A (v) = 0. . 

The criticality condition in this case can le reduced to 
the fajniliar ' end-point result ' [ 5 ] . 

h = b , . ^ |v j - z 

c end-point 2 ' o' o 

The corresponding solution is exactly of the same form as 
diffusion theory except the use of transport theory Vq 
place of asymptotic decay constant. In the first-order 
approximation, A'(v) in equation (1.5) was set equal to the 
inhomogeneous term A^ (v) and the corresponding solution was 
obtained after fairly lengthy calculations. Although, in 
principle, the scheme is applicable for higher order approxi- 
mations, no such attempt was made by Mitsis, because of the 
complicated nature of functions involved. It is to be 
noted, however, tha.t the Ileumann series converges rapidly 
except for systems with dimensions less than one mean free 
path i.e, for large c. Thus, for critical thickness calcu- 
lation the first order or even the zero-th order approximation 
is fairly accurate for all practical purposes. In studying 
the slab albedo problem, McCormick and Mendelson [17] 
followed an approach similar to that of Mitsis and constructed 
a slightly different iterative scheme for the solution of the 
coupled Fredholm equations based on Neumann series expansion. 
The Albedo problem transforms to the criticality problem in 
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ulie absence of incident neutrons a,nd tlie results obtained 
by HcCormik and Mendelson are same as Mitsis’s. The 
conser^vative convergence criterion for Neumann series that 
the modulus of eigenvalue bo greater than the spectral 
ladius of the opeiator, puts a restriction on the validity 
of solution for the Predholm equations. Thus the inte- 
gral equations are not solvable by Neumann series expansion 
for a certain range of parameters although the solution 
actually exists in those cases. Both Leonard and Mullikin 
[18] and Bowden et al [19] apply a similar approach to 
construct Neumann series solutions for their integral 
equations, Erdmann [20] gives an elaborate discussion on 
different approximation schemes for Neumann series solution 
of the coupled integral equations that arise in the case 
of neutron transport problems in finite media. In none of 
the works mentioned sofar, the integral equations were 
actually solved - zero-th order or at best first order 
approximations could onljr be constructed. 

Among the approximate methods, the spherical 
harmonics approximation is simple and popular as an 
analytical tool, whereas the discrete angle approximation 
is more suited for numerical work, A good accoimt of 
these .methods are compiled in reference [6], We concen- 
trate on Pjj- approximation only as it is of relevance to 
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the present ■work. The detailed analysis of Pjj-approximation 
can be found in two classical pa.pers by I%rsha,k [21] a-nd 
I'dark [22], and the t'frJ’o widely-used free surface bounda.ry 
conditions in this approximate method go by their name. 

The basic idea of the method is to expand the neutron 
angular flux in a series of Legendre polynomial and replaces 
the transport equa,tion by an infinite set of coupled 
differential equations in terms of the moments of angular 
neutron density. In approximation, the infinite set 
is truncated to a finite set of (L+l) equations by imposing 
the somewhat arbitrary condition that the derivative of 
(F+l)-th moment of the angular flux vanishes. Because of 
this truncation, the continuous spectrum v e (-1, 1) is 
replaced by (N - l) discrete values and thus the eigen- 
value spectrum is purely discrete. A major short coming of 
conventional approxima,tion is that for small values of 
c and 5r not too large, all the (N + 1) eigenvalues . lie in 
C-1, 1) while it is known that the largest ) v.j corres- 
ponding to the as 3 nnptotic flux is greater than 1, The 
error involved in this approximation is mainly due to the 
truncation of infinite number of differential equations to 
a finite set and also due to the fact that exact boimdary 
condition cannot be satisfied in such an approximation. 

Ever since the approximation was proposed for neutron 
transport ealculation, much effort has been made towards 
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its improvement. Tlie first such attempt was on the or 
classical diffusion appro xima.tion. This improved theory, 
C8,lled the asymptotic diffusion theory, is described by 
Case, de Hoffmann and Placzek [23] and by Davison [4]. 

In this theory, the diffusion coefficient of P^ approsima- 

p 

tion is replaced by (l-c) v /E on an adhoc basis. Although 
the exponential growth or decay of neutrons flux is 
accurately described by this theory, the stringent boundary 
condition of continuity of scalar flux at an interface 
lea.ds to wrong flux amplitude deep within the medium. 

Pomraning and Clark [24] developed a slightly different 
asymptotic diffusion theory based on variational methods 
that gives exact transport growth or decay and also 
predicts the flux amplitude correctly. This was done by 
allowing a flux discontinuity at an interface. Selengut [25] 
also used variational methods to introduce the as 3 niiptotic 
diffusions theory where it is foimd that both current and 
scalar flux are discontinuous at an interface. Pomraning [26] 
later developed another a.symptotic diffusion theory to 
incorporate the advantages of both iiie earlier methods of 
Pomraning and Clark [24] and Selengut [25], In his 
asymptotic theory, the scalar flux is allowed to be disconti- 
nuous by introduction of two parameters p(c) and r(c). 

These parameters are eva.luated by normalizing the formalism 
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to two standard reference half-space prohlem solutions, A 
supplementary separate diffusion-like theory is proposed 
for tia,nsicnt flux. Unlike the asymptotic flux, the total 
(asymptotic plus transient) scalar flux is continuous 
a.cross a hoimdary. This is achieved hy introducing another 
parameter a(c) in the tra.nsient flux theory which is deter- 
mined from the condition that the total flux gives exact 
transport-theory leakage from a constant source Milne 
problem. This excellent theory proposed by Pomraning has, 
however, certain limitations. Although boundary conditions 
have been modified through introduction of different 
parameters, a generalized modified boimdary condition has not 
been prescribed. This wa.s not possible because the forma- 
lism did not go beyond a diffusion-like approximation. Also, 
the theory is based on replacement of vacuum by purely 
absorbing medium (i.c. Hark boundary conditions) - the more 
versatile Marshak boundary conditions were not used. The 
theory is only concerned with scalar flux - thus improve- 
ments on various integral parameters (angle-integrated) 
were only tested. Some inconsistency is also observed in 
evaluation of r (c) for c > 1, and a(c) cannot be calculated 
for c > 1. Calamo [27] prescribed alternate sets of para- 
meters for use in Pomraning's diffusion— transport equation 
to give flux transients at an interface and near source 
accurately. Doyas and Koponen [28] extended the asymptotic 



neutron diffusion theory to linear anisotropic case. They 
suggested and evalimted anisotropic transport-theory 
corrections to 8.symptotic diffusion theory. The improve- 
ment of general Pj^-approximation along with the work on 
as 5 niiptotic diffusion theory was also in progress. 

Conkie [29] gave several iterative schemes to improve 
the conventional Pj^-approsimation, In his analysis, itera- 
tion procedure started with ^the Pj^-approximation -ra-lue of 
the angle-integrated flux. The idea was that the total 
flux is determined more accurately than the hi^er angular 
moments, since it is an integral over the approximate 
angular distribution and presumably the error in the angu- 
lar distribution is oscilatory because of orthogonal 
function expansion. It was also possible to satisfy the 
exact boundary condition of the problem. In his later 
paper, 'Pomraning [50] presented an approximation to the 
transport equation which yields the exact transport theory 
asymptotic behaviour in all orders for any geometry and 
is accurate to any arbitrary degree. The approximation, 
very similar to the usual spherical harmonics method is 
derived by introducing a new truncation scheme into the 
infinite set of spherical harmonics equations. The trunca- 
tion method consists of assuming that the higher spherical 
harmonics components, equated to zero in the Pjj-method, can 
be related to lower components by assuming the angular 
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distribution to be in an asymptotic distribution. The 
resulting approximation, called as approximation, 
is not more complex than the usual spherical harmonics 
method for computational purposes. The (N+l)-th spherical 
harmonics equation in slab geometry is given as 


d , (x) d (j) (x) 

( 1 + 1 ) + (2n+l) (x) + N = 0 


(1-4) 


where . 


r 

<!> (x) = 2 % dp (x, p) 


(1.5) 


If the angular distribution ^ (x, p) in equation (1.5) 

is assumed to be in its asymptotic form as 


4 > (x, p) 


1 + p/ V 1 ~ p/ V 


( 1 . 6 ) 


where A and B are arbitrary constants. By use of equations 
(1.5) and (1.6) a relation between i^^^ 

derived as 


^H+l = “if 


(1.7) 


Where, 


% - 


I dp P 
J-1 
1 


1+1 (p) 
dp Pjj_^ (p) 1(2 .+p7vT 
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in equation (1.4) 
Introducing a parameter y-g 


% 


, /IT+l. -- 


is now replaced by relation (1.7). 
as, 


this process is equivalent to setting equal to 

zero and replacing <{> "by Yjj equation 

(1.4). Copic [51] extended Pomraning's work to give an 
explicit expression for the ratio of (lir+l)-th to (K'-l)~th 
angular moments i.e. a^j. following Davison [4], he also 
showed that the attenuation distances of the remaining 
solutions (i.e, the transient component) in the finite 
order approximation are all shorter than the diffusion 
length as given by the asymptotic eigenvalue. The eigen- 
values in the Ajjj— approximation are given by a certain 
determinental condition just like the case of usual spheri- 
cal harmonics method. Although ihe largest eigenvalue 
gives the exact asymptotic behavior, it has not been shown 
in Pomraning's theory that the remaining ^ (H-l) pair 
of eigenvalues for odd order approximation are confined 
in the interval (-1, 1). The behaviour of A^j-approximation 
as IT “ was not also investigated.' An analytic transient 
flux expression for use with asymptotic diffusion theory is 
prescribed by Beghian [52] in his treatment of monoenergetic 
integral transport equation. The exponential transient 
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representation is replaced by exponential integrals' to 
yield significantly improved results over diffusion theory. 
Summarising, it can be said that the improvement on conven- 
tional Pjj approximation proceeded along t'wo parallel paths: 
one was the development of asymptotic diffusion theory by 
replacing the diffusion parameters like decay constant, 
diffusion coefficient by their transport theory values. This 
was mostly achieved by either ad-hoc intuitive basis or 
variational formulation. The other was on development of 
improved general P|j— theory with' transient components 
included. Unfortunately, the development in this direction 
lacks an unified and consistent approach. In many cases, 
systematic construction of higher order approximation is 
not possible. Most of the theories are concerned with 
scalar flux except perhaps Pomraning’s theory which deals 
with angular flux. However, the theorj?" has its limitations 
as pointed out earlier. 

Outlin e an d O bj ect ive of T hesis : 

It has been outlined in the ' previous section that 
the present work on one-speed neutron transport theory is 
centered around the spherical harmonics method and Case’s 
singular eigenfunction expansion method. It is true that 
the Pjj-approximation cannot be a competitor to the versa- 
tile Sjj-method for practical reactor calculations , but the 
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former method has some plus points in some cases. The 
heauty of Pj^-approximation is that it can provide the 
functional behaviour of angular neutron density while the 

Sjyj method deals with discrete numbers only. Due to this 

% 

factor, the old spherical harmonics method has regained 
its lost impetus in recent times- as for example, the 
prediction of functional behaviour? of neutron angular 
density deep within a medium by Kin [33] in his disserta- 
tion. The improvement of conventional Pj^-approximation 
is necessiated by faster convergence of solution. Prom 
what has been sketched in the previous section, it can 
be concluded that while the state of art has proceeded 
remarkably in scalar flux improvement by replacement 
of diffusion parameters by their transport theory counter- 
parts and by allowing some discontinuity in flux or 
current or in both at an interface, the general improve- 
ment of Pjj-approximation through better representation of 
angular flux by modification of usual boundary conditions 
has not received proper attention. The present w ork is 
an attempt towards a better angular flux representation. 
The main objective in this work is i) to demonstrate that 
the use of all the roots of Pjj- approximation is not ma, the- 
matically justified and hence, by implication, it is the 
transport theoretic which should be incoi^o rated in 
place of asymptotic decay constant of Pjj--app3:‘oximation, 
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ii) to generalize the conventional Marshal free-surface 
boundary condition in order to get the fullest advantage 
of the use of exa.ct v^. As it will be seen in the 
3.118,1 ysis 5 the first part is achieved by suitable construc- 
tion of gj^(v) functions as they appear in theory both 
in the ranges of v i ("Ij 1) and v e (-1, l). Modification 
of boundary conditions is realized by introduction of a 
parameter cCj^Cc) in the first equation of the set of 
I'larshak boundary conditions. ctjj(c) is evaluated by forcing 
the asymptotic flux for the source- free Milne problem to be 
identically same as exact transport theory asymptotic flux, 
The Milne problem is chosen as a reference problem as it 
leads to the natural freo-surface boundary conditions. The 
same can be used for any other transport problem involving 
a free surface. The solution is constructed with the 
asymptotic part being of the same nature as exact solution 
and the tra,nsient part terminated as in conventional P^-appro- 
ximation. Since the transient is in the framework of 
Pjj- approximation, unlike previous attempts, construction 
of any order of approximation is possible in the present 
approach. Interface boundary conditions, which are more 
useful for practical reactor calculations, are also modified 
in the proposed TPjj theory. In a two-media Milne problem 
the modification is achieved by demanding an exact asympto- 
tic behaviour in the right hand side when the adjacent 
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medium is replaced by a purely absorbing one. It is 
expected in general that through exact angular represen- 
tation of fluX; atleast for the asymptotic part, the 
convergence of solution can be achieved in a reasonably 
lower order of approximation. However, due to inadequate 
transient flux representation, the TPjj-approximation does 
not necessarily improve the Pjj result in those cases where 
transient is overwhelmingly predominant.. This calls for 
a better transient flux representation and it is felt tha.t 
this may not be possible within the framework of Pjj-approxi- 
mations. Finally, it is to be concluded that the different 
inferences drawn in the present analysis, like replacement 
of different Pj^-approximation parameters by their transport- 

i 

theory counterparts, are not just adhoc prescriptions but 
are justified by rigorous mathematical formulation. 

r 

As has been remarked in the last section, all i 

i 

standard one-speed transport theory problems have been 
ana-lyzed by Case's exact theory in great detail. Thus, | 

from ca.lculational point of view, it is not worthwhile to I 

compete with the existing methods of solving these problems. | 
It is only some general mathematical problems associated | 

with Case's method that need attention. In neutron trans- j 

port theory, eigenvalue problems occur repeatedly in | 

different contexts where the existence of dominant eigenvalue | 

■ ■ ■ ' .. ■ 

^ 
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is of great importance. Proper mathematical study of such 
a problem was pioneered by Lehner and Wing [54] in their 
analysis of one-speed time dependent transport equation 
in a finite slab. The equation is transformed into an 
eigenvalue equation of the form 

Y r ^ 

<|)(x) = E.(^lx-x’I)<t)(x')dx' ' (1.8) 

-a 

Equation (1,8) can be viewed as an inverse eigenvalue 
equation in the sense tha,t the eigenvalue y/2 is known 
and the parameter X (spectrum of one-speed operator) which 
has to be determined is in the kernel of the integral 
equa.tion. The task would have been easier if A were 
fixed and y/^ were .to bo determined! By means of classi- 
cal approach, Leliner and Wing were able to prove the 
following result - the Boltzmann operator decomposes the 
X - plane into i) point spectrum - a finite non-empty set 
IjHng on X > 0 ii) continuous spectrum - the real line 

Re ( x.) _< 0, iii) residual spectrum - empty, iv) resolvent 
set Re (x) > 0 deleted by the point spectrum. The exis- 
tence of a dominant eigenvalue was established by estab- 
lishing a proper ordering among the discrete eigenvalues 
of the point spectrum. In his study of time dependent 
slowing down equation (spa.ce independent) Ukai [35] gets 
an equation similar to equation (1,8) of Lehner and Wing as 
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P 4> (v) 



G-X v') ^ 


(1.9) 


where v is the velocity variable. Ukai proceeds by fixing 
^ arbitrarily and solves equation (1.9) to obtain a set 
of eigenvalues, p^. Each eigenvalue P contains ^ as a 
parameter so that p is denoted as p (^). Prom equation 
(1.9) it is evident tha.t the parameter ^ is given by the 
solution of the equation o(^) = 1. Punctional analytic 

studies on general Boltzmann operator A were also made 
by several authors [36-38], In this case, it was found that 

i) the whole half -plane Re (x) <. “ ^* = ^ (^» ’'^)] 

belongs to continuous spectrum of the Boltzmann operator, 

ii) there exists a number c^ such that the half-plane 
Re (x) > Cq lies in the resolvent set of A and iii) the 
points of the strip - X* < Re ( x) <_ c^ also belongs to the 
resolvent set p (A) with the exception of a set of isolated 
eigenvalues which constitute the point spectrum. In his 
analysis, Vidav [36] applied a theorem from functional 
analysis to establish the existence of a fundamental mode 
of decay if the set of ^ ^ Is not empty. The existence of 
corresponding non-negative unique eigenfunction was proved 
by application of a theorem by Krein and Rutmann [39]. 
Similar analysis in a more general way can be found in a 
paper by Larsen and Zweifel [37]. In Chapter III of the 
present work, the two coupled Predholm equations in terms of 
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the discrete and continuous coefficients arising in a finite 
slab criticality problem are coupled together to form a 
single integral equation. The problem is then posed as an 
inverse eigenvalue problem in which the eigenvalue is given 
and the operator ha,s to be determined. By application 
of functional analytic theorems in the r ange of analyticity 
of the operator-valued function |f<^(b) (b is the inverse 
eigenvalue parameter), the existence of critical thickness 
is proved. Analyticity criterion on IK(^) shows that the 
critical thickness b^ must be less than the end-point result. 
This fact is not reflected by Pj^-calculation. Estimates of 
spectral radius in its simplest form are used for numerical 
calculation of critical thickness. It was necessary to 
depart from the original problem slightly because of 
restriction of positivity criterion on the kernel of the 
integral equations. This is because the relevant theorems 
are valid for positive operators only. As calculations 
will reveal, the results obtained are fairly accurate in 
spite of dealing with the approximate problem. The choice 
of the criticality problem was made as this problem has 
not been seen in this li^t. Also, the present analysis is 
based on a. more broad-based mathematical theory. The 
numerical results in this work were obtained only to subs- 
tantiate the hypothesis adopted in the formulation. Ho 
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sophist ication in calculational tecliniq,ue has been attempted 
to improve the numerical results, because our main aim was 
not to calculate critical thickness (which has been done 
to great accuracies earlier) but to formulate and solve a 
class of inverse eigenvalue problems. Thus the present 
calculational method does not claim to be a competitor to 
the fast converging Neumann series solution for critical 
thickness. However, in those cases where the numerical 
values of inverse eigenvalue parameter are not otherwise 
easily available, the present method of calculation may be 
adopted. 



CKiPTBE II 


A TRANSPORT THEORETIC P^j-APPROXIIATION 

I ntro d uction 

In Chapter I, an outline of the present state of 
Pjj-approxiination is given. As has "been noticed, most of 
the works on improvement of conventional Pjjj.-approximation 
deal with scalar flux 8-nd they are diffusion-like approxi- 
ma.tions. In this chapter, we give an improved approximate 
scheme that can yield letter angular flux representation 
and can he easily adopted for higher order approximations 
heyond P^^ approximation. To realise this, the discretization 
procedure of eigenvalue spectrum in conventional Ppj-approxi- 
mation is modified. This is necessiated hy the mathematical 
analysis that follows in the development of the first part 
of the theory. Simultaneously in the second part of the 
proposed theory, the usual free-surface boundary conditions 
are modified to achieve better angular flux representation. 
The improved scheme in genera,l can be applied for any 
problem involving a free surface. The usual interface 
boundary condition of continuity of moments of angular 
neutron density is also suitably changed in the proposed 
analysis. Formal solution of the transport eq^uation is 
then constructed with the asymptotic part being of the same 
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mture as exact solution and the transient terminated a.s 
in conventional Pjj-approximation, The complete formalism 
is called transport theoretic Pj^-approximation and this 
approximation will he referred as TPj^ approximation. This 
chapter ends with numerical calculation of different 
standard one-speed neutron transport problems by TPj^ method. 

2.1 Analysis 

The one -speed transport equation in slab-geometry 
with isotropic scattering is 

M- + E u) = ^ r i>(x, p') dp' (2.1) 

9 X ^ J_3_ 

In Case's exact analysis,- the solution of equation (2,1) 
can be written when c 1, as 


'I' (x, p) 


-ExA 


^0+ ® 


° + "o- 

+J d, 


Ex/v 

® ° 'J’ 


( 2 . 2 ) 


where. 


<Po+ = 


C V 


0 

Vo ? p 


( 2 . 3 ) 


in equation (2.3) is a solution of 


In - 1 

2 v-1 ~ 


(2,4) 
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and, if V e (-1, l) 

6 (p.) = 2^ P + A (v) 6 (v - i.i) (2.5) 

In equation (2.5), P denotes the Cauchy principal value. 

The main difference between equations (2.3) and (2.5), 
besides the principal value, is the appearance of the 
extra-factor involving the arbitrary function A (v). This 
fact will prove to be important in the analysis to follow 
on the spherical harmonics approximation. 

In the spherical harmonics method ip(x, p.) is written, as 

<i)(x, p) = J (--fci) '/'^(x) P^(p) (2.6) 

n=o 

Assuming , tlie above becomes 

*(x, |i) = J (^) g^(e P^(n) (2.7) 

n=o 

In exact treatment the infinite sum in equation (2.7) is 
denoted by <{>y (p) i.e., 

♦v (l-i) = f (--2-) Sa('’> (2.8) 

n=o 

The coefficients g^(v) in the expansion of (i>^(p) can be 
obtained simply as 

Sii(v) = f <^y(p) «3p (2.9) 
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By use of equations (2,3) and (2.5), we get from eqn.(2,9) 

g^_(v) = cvQ^(v), V i (-1, 1) (2.10) 

= CvPQ,^ (v) + X(v) Pji(v), ve (-1, 1) (2.11) 

The functions (v) are the Legendre polynomials of second 
kind and are defined as 

,1 P (ji) 

Q3^(v) =f \ d|a, V i (-1, 1) (2.12) 

^-1 (v~ ti) 

and , 

P Qj,(v) = I P J ^ ^ ^ (2.13) 

When V £ (-1, 1), principal value integration is used .Because 
of singularity at iJ, = . For this range of '* » functions 

will he denoted as PQj^(v). A Brief discussion on these 
functions is pertinent to the present analysis. The func- 
tions Qj^(v) like Pj^(v) a.re oscillatory in the interval 
(-3., 1) and possess zeroes there, But like P^ (v), they 
are monotonic functions outside this interval. For example, 
when n = 0 

QqC'') = In “J r V ^ (-1, 1) 

and , 

Qj^(v) = I In lyy , V £ (-1, 1) 
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It is seen tha,t Qq(^) po.ssesses zero at v = 0 for 
V e (-1, 1), but does not have a zero for v ^ (-1, 1). 

The Q^^(v) functions satisfy the same recurrence relation 
as the v) , 

(n-i-1) = (2n+l) vQ^^(v) - n 1 

(2.14) 

With the knowledge of Qq(^) and Q 2 _(v)» for n >1 

can he evaluated from equation (2,14). Qq(v) can be 
obtained from equation (2.12) or (2.15) as 

Ql(v) = vQ^(v) - 1 V ^ (-1, 1) 

VP ;Q^(v) - 1 V e (-1, 1) 

Qj^(v) can also be expressed in terms of 5j^(v) by the follow- 
ing continued fraction expression 

, 114 (n-1)^ 

rjvr = - VI 3U 5— ••• - Tz'n'-'iiv 

Another interesting property of Qj^(’^) is that Qj 3 ^('') 0 as 

n ~ for all v , 

r ^ 

The normalization condition ; *!*u(^'-) d|j, = 1 for 

-'‘-I 

all V implies g^i^) = 1. For v y (-1, 1), this gives 
Case’s tianscendental equation (2.4) the solution 
of which are the pair of asymptotic decay constants + 
while for e (- 1, l), this determines the foimi of v 
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dependence of ^ (v)“ 

X (v) = 1 - c V P (v) 

With this value of X(v), the expression for (v ) in the 
range v e (-1, 1) is that given in Davison [4]. Plowever, 
the present analysis shows that this expression for 
is valid only when ve (-1, 1). The conventional P|;j-“apprO” 
ximation involves approximation of equation (2,8) by a 
finite sum i»e., 

g]jf+l (v) = 0 (2.15) 


The equation (2,15) forms the basis of Pjj--approximation, 

The general solution in such an approximation is constructed as 


4(x, p) = I 


3=roo bS 


a . ^ " e 

j 


(M) .-“/"j f (S-l) g^(v )P^(^) 

n=:0 

(2.16) 


where the entire range of v is suitably discretised, 


It is obvious that gj^(v ) of eqn, (2,8) cannot be zero for 
any finite n = IT because the Legendre polsmomials are an 
infinite set of orthogonal polynomials. (Eqns. (2,10) and 
(2,11) ofcourse satisfy the condition = o as n -»• “>). 

Any attempt to set S^j+qCh') = 0 can, therefore, result in 

possible inconsistencies. Thus, Davison [4] obtains the 
largest root (aS 3 nnptotic) ' less than 1 for small values of 
c. Looking from equations (2,10) and (2.11), the above 
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follows in a straightforwai'd manner. When v is outside the 
interval (-1, l), 8rx('^) given by equation (2,10), The 
properties of show tliat it is impossible to satisfy 

condition (2.15) when i (-Ij 1). Thus the roots outside 
(-1, l) (i.e, the asymptotic roots) can never be found 
from the condition (2,15) and must always be evaluated from 
the normalization condition 

§0 = ^^o '^o =1 ^2.17) 

In the range of v s (-1, 1), it can be shown with the proper- 
ties of Q^(v) functions that the expression for 
equation (2,11) reduces to a polynomial of (lT+l)-th order. 

This poljmomial must have (N+l)/2 pair of zeroes and thus 
the condition (2,15) can be satisfied in the range of 
^ s (-1, 1). However, since equation (2.15) in this case 
has the required (N+l) roots, two of them must be redundant 
as they have been already obtained from equation (2.17). It 
is known that equation (2.15) can admit largest pair of 
roots both gre-ater or less than 1 [4]. When greater than 1, 
this pair is straightaway rejected as equation (2,11) is 
valid only for v e (-1, 1). However, the inclusion of the 
largest root when it is less than 1 is not also possible 
as this would lead to consideration of two extra roots. We, 
therefore, conclude that, irrespective of the value of 
largest roots', only the smallest (N-l) roots of equation (2,15) 
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(where eqn. (2,3.1) is used for ^(v)) should he used fo 
describe the transient flux while the asymptotic flux has 
to be constructed from the roots obtained from eqn. (2,17). 
In his analysis, Hendry [40] ne^qlects the roots greater 
than 1 but retains them as two extra 'basis functions' when 
they are in the transient range. This is, however, not 
true as the two additional roots do not generate extra, 
basis functions; the orthonormal basis functions in the 
present case are the Legendre polynomials and not the 
spatial modes exp (-x/v). The present analysis also 
establishes an exact relationship between (H+l)--th moment 
to (N~l)~th moment of angular flux, a quantity used by 
Pomraning [30] in his A-^x approximation theory. The parame- 
ter ajj, as defined by Pomraning is 




H+1 

N-1 


(v) 

'(7) 


Prom eqn. (2,8). it follows, 

%+l(v) 


ajj _ 




%+ l(^) 

%_l(^) 


by definition. Thus, the replacement of by . 

‘J'jj ^(x) in F-th eqn, of Pjj.-approximation will obviously 
lead to correct asymptotic behaviour of neutron flux. 
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In Appendix, functions are constructed for 

anisotropic scattering. 

free Su rface Boundary Cond iti on; 

The replacement of va.lue of the asjnnptotic decay 
constant hy its exact value modifies the Pj^-solution to 
(e = 1), 


(3C, \X) = 


(N) 

* 0 + 




N 


° I g,(v,) 


n=o 


a») P,(.) 

n=o 


I 

j=roots 


W I 


n :=0 


I P,(u) 


where v are (K-l) roots corresponding to the transient 

J - 

flux and are obtained from solutions of eq,ns,(2.1l) and (2.15) 
However, considering the poor convergence of the legendre 
Polynomial series in n,[4l] we wish to go further and thus 
retain the infinite sums for the asymptotic roots and 
terminate the series for the transient part only. Hence 
the general solution changes as 

♦ ,) =a(f -aW e"''''" 


I e I g„(v .)P (n) 


j=roots 


n=o 


n' j' n' 

( 2 . 18 ) 
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The coefficients in eqn. (2.18) are to he evaluated from 
the boundary conditions. Equation (2,18) shows that asymp 
to tic part of the flux in this case is of exact nature, 
and can indeed be made equal to its exact value by 
suitable modification of boundary conditions. But in 
order to do so, the boundary condition has to be changed 
for each problem to ensure exact asymptotic flux in each 
case. This idea is not pra.ctically applicable. To gene- 
ralize concept, we wish to get the exact extrapolated end 
point and modify the free surface I-Iarshak boundary condi- 
tion, Since the source-free Milne problem defines the 
extrapolation end-point , where the asymptotic scalar flux 
vanishes, the modification of the I'larshak boundary condi- 
tions at the free surface will make the asymptotic 
coefficient equal to its exact value for the 

Milne problem. The modified form of tiarshak boundary 
condition is then applied for any other problem involving 
a freo-surface. While for source-free Milne problem the 
asymptotic flux will be exact by definition, the same will 
not be true for other problems but will d ofinitely be an 
improvement. This conclusion will be reflected by 
numerical calculations of different standard one-speed 
problems given at the end of the chapter. 
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2 • 3 So ur c e-Fr e _e_ Problem and Mo d i fl ed Bo undary 

Conditions 


The exs.ct solution of the source-free Milne 
problem is 


-s:/ V x/v 


<!> 




A(v) 

n 


-x/v / \ I 
e 6 (p) dv 


(2U9) 


where , 


0 + 


I(v„) 


•2z /v 
o' o 


( 2 . 20 ) 


where X(v^) is the well-known Z - function of Case's theory 
and is defined by the eqn, [42] 


Z(z) = 


2.^ . f __1L1 - 

2 l-G ii (v^2 _ ^,2) x(p') ii'+z 


z^(c) is the exact extrapolated end point and is 
tabulated by case, de Hoffmann and Placzek [25]. 

The associated free-surface boundary condition is 

A(v) <^y(p) dv , Ji > 0 

( 2 . 21 ) 
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In the present case, the exact solution as given by 
eq.n, (2,19) changes to 


’i'(x, \i) 


-x/v, 




4>o+(lJ-) + 


x/v 


0o.(p-) 


I’ 


-j - ^ > i -)%cv4)p„c^) 

3 =roots J n:;=o ^ ^ 

( 2 , 22 ) 


The associated free-surface condition is 


^ ° ^ 3 =roots n=o j ix 


( 2 . 23 ) 


where, prime in the summation indicates that only positive 
roots are considered. 


The condition (2,21) is obtained from the exact free— 
sui-face conditions 4'(x, lOt^ji-o ~ ?. ^i^ce ^^^^(p) 

and (sj.) form a complete set in case of half-space 
problem, a^^ and A(v) of oqn. (2,21) can be exactly evaluated 
from half • range orthogonality relations. Because of 
trimcation in the transient part, eqn. (2,23) cannot, 
however, beexactly satisfied. In the conventional Marshak 
boundary condition, moments (n = 1,3,.. , S') of ^(0, ji.) 
are only made equal to zero and these equations are used 
to evaluate the coefficients. Since of eqn, (2,23) 



31 


are infinite sei-ies in Legendre polynomials ( eqn. ( 2,8) ) , 
the Legendre polsnaomials are still the basis functions j-n 
the. present case. Therefore, in the modified version of 
the free surface boundary condition the moments of 

ii) are to be used. However, in the present forma- 
lism the number of unknowns are reduced to (H-l) as the 
asymptotic coefficient is forced to its exact value. Thus, 
it is necessary to introduce a parameter a^^Cc) in the 
first of the set of I'fe.rshak boundary conditions so that 
a^^ is indeed given by eqn.(2,20). a|^j(c) will be evaluated 
from the condition of exactness of asymptotic coefficient. 
The detailed analysis is given below. The conventional 
]>Iarshak boundary condition is 


( P^(p) ^(o, p) dn =0 n = 1,3,5... H 

Jo 

The modified free-surface condition becomes, 
f (p + ocjj) ’l'(0, ij.) dp = 0 

Jq 


(2.24) 


(2.25a) 


C Pj,(li) ’>'(0, p) dp = 0 n=3,5...H (2.25b) 

: 0 

i.e. instea,d of taking {P„(p) ^ „>moments, we 

take here {(p+a), n=3 5 . ^ “foments of t(.0, p). 

The use of (2, 25) in eqn. (2,23) ^delds, 
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C V. 


'1 ( [ J- -f OCj^) 


\ d^i - I - a 

-^O (v^+ li) °+ 


(\1 + a^) 
o (vq - !-i) 


d^i 


and 


' fi.T'i A 

+ I J^(2n+1) ° 


(2.26a) 




HI 

c ) — -^ — — (ip^ _i- 0 ^^ ^ 


■1 




JT JL 

I (2n+i) gj,(v.)f p (n) P (u) an = 
3 n=o J J o 


o+ " "o j '7 ”‘7 
-^o 


0 


m = 3,5,...N (2.26 Td) 

Equations (2,26a) and (2,26b) are (“■^) eo^uatior^in the 
unhnown and (^^) coefficients It is to be 

O 

noted that aQ_j_ in this case is known and is given by 
eqn. (2.20). thus obtained is used for all problems 
involving ha,lf-space. The formalism is now complete and 
the present scheme willj henceforth,; be called as TPu 
approximation. In the next section, the evaluation of 
for TP 2 _ and TP^ are presented. 

2.4 a|j(c) for TP^ and TP^ Approximations 
In the TP^ approximation, we have 
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i (|j. + a, ) 

(Jn + a . 

Jo (» + n) 


J o 


1 (p. + a, ) 


(v^ - ) In 


X. e. a 


0"1” 


(^o 


V, 


o 


By use of eqn* (2*20) ^ we get. 


or, 


-e 


o o 


(^o ~ 


-z /v 
5 O O 

■■rr;7v7 

o o 


-H ) In 


0 

- p) 

1 


1 

- 1 


- 1 

1 


V 

+ 1 




^o 


^'0 

^0 

1 


+ 1 

0 


0 


V 


0 


dp, = 0 


- 1 


- 1 


- 1 


^ 1 


z /v — z /v 

0 o 0 0 

e - e 

z A "" -2 /v 

o' o 0^0 

e + e 


V +1 

xT 

o 


V 

o 

o 


a-. In + V In 

1 v^-l 0 


jiOr c < 1, the left hand side is tanh (^qAq) and for 
c >1, using Vq = ^®aomes tan (z^/k^). Also, 

hy definition the linear extrapolation length X (c) 
is given as. 
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x(c) = tanli c < 1 


tan c > 1 

o o o 


Thus, when c <1, we get 


A(c) = 

i 


2v ^(l^c) + a-, c ^ ^ In - -|-— 
® J- o ^ 0-^1 


2cx-j + c In — ^ S 
1 O V ^-1 

0 


and when c > 1< 


(c) = 


2 

P ^ k "fl 

2k^ (l"c) - a3_ c In - 


P + 1 

2cc-, c k In , - 

1 o k^ 

0 


a-, (c) is obtained from (2.27) and (2.28) as 


a. (c) 


2(l^c) - c^{o) In — p’-^“ 

"o-^- 


C< ] 


2^(c) - cvf In 

2 1 

V «-l 

o 


and. 


a^Cc) = k; 


2 

k^ + 1 

2(c-l) - cX(c) In -p- ■ 

k"^ 

p___, ^ , 

, . 2 k^ +1 

2A(c) - ck^ In ■ 


(2,27) 


(2.28) 


(2.29) 


1 (2.50) 
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In TP approximation, the boundary conditions in our 

3 

•-formalism are. 


fl ( ) 

c V. 1 du -F a . cv. 


''<=1 


o (v^+p.) 


0 + o 


^1 ( ) 

O (^ Q - P -) 


dp. 


1 


+ I (2n+l) g (v^)i (p+a^) <3|i = 

n=o Jo 


= 0 
(2.31) 


and. 


fl P^Cp-) 

°Jo T^) 


/I P-T ( [J. ) 

cv_. — « d|i -I- — -»d|i 


o- 


° J 0 


+ Ap^ I (2n+l) g,,(^T); P,,(p) P,(p) dp = 0 (2, 


n=o 


1.^1 n 
^ n 


fl ?„(],>) 

To evaluate the integrals of the tjrpe < — — ' - dp , follow- 

Jo (^Q-'P) 

ing results are used [43] 


(n) 

= J 


dp 

1 

■+ 

Jo 

( Vo + H ) 

n 

(n) 

H O 

ti 

n 

.....ih 

(Vq-u) 

dp 

_ _ 1 
“ n 

(0) 

•+ 

= In 

Vi 

V ^ 

and 

po) 




(l-nv^ I 


= In 


(n-l) 


o 


V _1 • 

o 


32) 


with, 
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Thus - 


A P^Ci^) 

(^o ± 1J-) 


5 5. T 

2 + “2 ■^+ 


Carrying out the integrations in eqns, (2,31) and (2.32), 
these eq.ua.tions simplify to 


= '-0 + <^3 1^°’) + s+ “'■o + “3 


+ Ap^ { {-j + 6\('']_) + I S2(''i)) 


+ (1 + “ - Q (’^p)) ^ 


(2.33) 


-o (I I?’ - i -0 (| - i ^1"’) 


+ Ap^ {- g + I ^ ° 


(2.34) 


Eliminating An from eqns, (2.33) and (2.34), the expression 

-2z /v 

for a^(c) is obtained as, with a^^ = -e 


a3.(c) = 


%2 ^'^l^ PpCc, ^q) + ^23*''^!^ ^2 ''o^ 

^23 =^3^°» '^o^ 1^31 ^'^l^ ^4 ''o.^ 

(2.35) 


c < 1 


“^(c) 


^12 ^'^l^ ^o^ ^23 ^’^1^ ^o^ 


( 2 . 36 ) 


, c > 1 
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where, 


and 


p, 

i^l 2 


('^l) - 4 ig" S2(''2.) 


P23 (^ ) = - 


16 16 2 ^3 ^'1 


-i- ) + — g (v ) (2^71 


hi “ 


1 , 3 


^ 4 - 16 h 


V 


f^Co, \) 


o 

2 


X (c) [3 In “2® 




5(v In — — 

O 2 n 

Vq -1 


-v„ { 3 (| - 1) - 5 (i -D- |]) 


^ 2 ^^» Vq) = 


x(c) in -2'-^- - 2V^ (| - 1) 
o 


^3(0, v^) 


V +1 

V In ^A- - - X(c) In -2--.- 

° vj -1 ^o"’^ 

o 


''o> 


v„ <3(i-l) -5 [v^ (i -1) 


2 ,1 

c 

2 


> 

3-' 


- X (c) [3 In -2^- 

^ r -1 

O 

V ^ 

- 5 (v^ In -1)] 

o 


.-2 


fl(c, k^) 


k? +1 


“#• ^(c) [- 3 In —2“ 
^ k 

0 


^ kf +1 

^ 5 (kj in ^-^2 

o 


- 1 )] 
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f 2 (c, ICq) = 
f 3 ( c ; ICq ) = 

f 4 ( c , = 

The procedure for evaluation, of ocjj(c) in the limiting cases 
c = 1 and c = 0 are demonstrated below. 

i) a 2 _(c) and a^(c) for the case c = 1 

For 0 = 1 ;, the two asymptotic roots coalese at 
-i- = ” and the two linearly independent eigenmodes are 

la ) =1 

'p2(:zy li) = ^ (-- !-!•) 

For source-free Milne problem, the general solution for 
c = 1 can be written as. 

♦(x, ij) = --f + i ^ e ^ ^ 

X f (-^a+l) g^(vj) ( 2 . 38 ) 


kf (-3(1 - i) - 5[M (1- i)- i]> 




+1 

0 


2 !^o ( 1 - 


2 + 2 - 
1 O 

k^ In — > 5 — 

• ° k ^ 

o 


2 X 


(c) 


■>^C ‘-2(1 - ?) - 5 [k^d- “) - j] ' 


k 

2 


k^+i 


2 - ^(c) [3 In - P -2 + 5 (k^ In ° 


k^-i -1 


k ‘ 


k 


- 1 )] 
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It is to be noted that z^, the exact extrapolation length, 
used in eqn, (2,38) ensures exact asymptotic flux for the 
source-free Mine problem. The free-surface boundary 
conditions are obtained from the conditions (2,25). Thus, 
in TP 2 _ approximation, the boundary condition is 




dia 


0 


or, a 3 _(l) 


3 Ml) - 2 

3(2z^-l) 3(2Ml)-l) 


(2.39) 


For c = 1, Zq( 1) = ^(1) = 0.7104 and a(l) is then evalimted 

from eqn. (2,39). 

For TP^-approximation the boundary conditions are 


■1 z -n, (3) 2 o T 

(ij. + a^) { J (^l+i) g^(v^) P^(!-i)> dju 

(2.40a) 


and 


f F^(ix) { + Ap) f g(v^)Fjix)} d|i = 0 

0 n=o 

(2,40b) 

(3) 

Eliminating from eqn, (2.40), a^(l) is obtained as. 


3Pi2(vi) X(l) + 4523 ^i) r3X(l)-2] 
~ 5p33_(v3_) x(l) +12 p 23 ('’iK 2 x(i)-.i] 


(2.41) 
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It is to be noted that in this case, gl( ) ~ ® 

X(l) = 0,7104. 

ii) aQ^(c) and a^(c) for c = 0 

The limiting values of cC]^(c) and a^(c) as c -*■ 0 3 ,re 
evaluated by a limiting procedure ^ 


As c 0, 


-i- = 

yj 

O 


By use of this, we get; 
2 


V 


In 


0 


o 


= ^ In (1 

O 


= £ -- In 4 

c 


So, from eq.n. 

a^(0) 


(2.29), with X(o) = 1, 

“2c -I- c In 4 
2-c(— I'- “ In 4) 


a^(0) is obtained as, 


2 “ In 4 
“~i5r? 


( 2 . 42 ) 


Similarly, in TP^ approximation, a^(0) is obtained by a 
limiting procedure from eqn, (2.55). 

PlpC^f) (In 4 - -J) + “ In 4) - 

^ 2 . 3 ^^ 111 4 + ^'^l^ (in ^ *” ^ ) 

The numerical values of Cj^Cc) for various values of c 
in (0,2) as obtained from ecin. {2,27)-(2.50) and (2.3S), 
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(2,41)-(2,43) are given in Table 2.1. Due to extreme 
iii"jensitivit,y of on c when c is very small, it is 
recommended that a^jCc) he obtained by linear interpolation 
between 0.0 _< c _< 0,2 in both TP^^ and TP^ approximations. 
Though it is possible to evaluate for higher order 
approximations (IT > 5) in a similar manner, the algebraic 
complexity of the expressions increases prohibitively. 

To get around this difficulty, we postulate that a^Cc) 
varies linearly with l/D between H = 3 and H = When 
IT ocu is zero as the TP^^ approximation tends to the 

exact solution in that case. Because the P^^ approximation, 
unlike all other higher order approximations, does not 
have a transient component, a-|^(c) is left out in this 
interpolation. Hence aj>j(c) is obtained as, 

^IT ^ ^ ^ ~ It ^3 (2, 44 ) 

2,5 Appli cations ; 

The TPjj approximation is now applied for the 
following neutron transport problems; the emergent angular 
flux distribution and leakage for source-free Milne-problem, 
the scalar flux for source-free Milne problem in case of 
c =1, the leakage in constant source Milne problem, 
the critical thickness for a slab and the disadvantage 
factor for a fuel-moderator cell, j^ongst the problems 
considered, the source-free Milne problem d efines the 




Table 

2.1; Values of 

as a function of c. 

■ c 

TP3_ 


c 

TP^ 

TP^ 

0.0 

-0.442695 

-0.522112 • 

1.2 

-0.036402 

-0.058201 

0.2 

-0.436811 

-0.514121 ! 

i 

1.4 

0.023245 

-0,069446 

0.4 

-0.375026 

-0.483447 j 

t 

1.6 

0.076578 

-0.063707 

0.6 

-0.273996 

-0„. 295818 1 

1.8 

0.125768 

-0.044994 

0.8 

-0.182374 

-0.076254 ; 

j 

2.0 

0.171308 

-0.018549 

1.0 

-0.103929 

-0.036359 : 






paramo t ex'* a-,:j-(c) by exactness of asymptotic fluxj while 
otliGi- pi'ohlems, which do not lead to the definition of 
the parameter, are considered to check the usefulness of 
the present formulation. 

2.5.1 Emergent Angular Eistrihution i{) (0 ,ijl) for Source- free 
Milne Problems 

By putting x = 0 in eqn. (2,22), ip(0, p) is 

p) = (p) + 

+ I’aW f s (v )r (n) (2.45) 

j n=o 

For approximation, as is obvious, the emergent distri- 
bution will be due to the contribution of exact asymptotic 
flux. It is given by, 

4'(0, p) = a^^ (p) + <!>^_(p) (2,46) 


In the TP^ case, for a source-free Milne problem the asymptotic 
coefficient is, by definition, eoLual to the exact value 

(3) 

a^^ and the only unknown coefficient to be deteimined is , 

This can be evaluated from equation (2,26b) with K == 3. Thus 
we get, 


A /I 

CV \ '7 - — V dp + a , GV i 7 < dp 

O Jq (Vq+pT o+ 0 (vQ-li) 

3 1 

+ Ap> [ (211+1) Sn(vj_)j P3(|l) Pj^(u) dn = 0 
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By simplification, is given as 


P 23 -^p 


(5) 


C V 


IT (5Ip^ - 


+ (51 


( 3 ) 


-f 


3 lj^hl 


and the emergent distribution is 


ij^( 0 , p) = a^, <}>„, (p) + 


(2.47) 


(3) 


+ AV' I g P (^) (2.48) 

n=o 


The exact angular distribution is [5] 

2 


i|;( 0 , p) = C’ 


X(-vp 


(vo - |i^)X(u) 


X(-Vq) is determined by solving the following two equations, 

^ (1-c) -1 


and 


X(v^) X (-V ) = 

O v_ (v 


-o ('-o -IV 


X * O ♦ y 


x77J 


X(-vp 


= a 


o^'• 


a 


= { 


04 - 


(1-c) -1 2 

° > 


^-7 ;y(;T.i) 


0 0 

In the conventional approximation, the emergent distri- 
bution is given by, 
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'('(O, ix) = I' i (™tl) g (v ) P (ij.) 

J j^_Q ■- n J ^ 


3 


K 


,(N) 


+ I (2n+l.) g^(-vy^O Pj,(|a) 

r.=o 
f IT ) 

where v ^ is Pj^-approximation asymptotic root* The 


coefficients q.\ 

vl 


(K) 


are determined from the usual Ifeirshak 
boundary conditions described by eqn.(2.24). Por approxi- 
mation, this gives, 

4I /-] ') 1 -t: (1) 

! ' H (^ + I gi (% )) <3|i 

o 

■ 1 


+ 


( (^ - I Ip. = 0 


or. 


(1) a 

“0 H 


a; ' (/I ■'■ 2 Si ) ) 




i.e. 


a 


( 1 ) 


1 

2, 

1 


( 1 ) 


2 o ^ 


Cl) 


^ + si(»ro 

/-j \ 

is the positive root of g2(’^) = 0 , with §2^’^^ given 
by ( 2.11) . Thus , 

^(1) ^ 1 

mr-cT 

By use of gp('^) = '^(l“c)s ^0+^ then becomes 


(1) l_ 


1 

2 


J a-c. 

3 


1 + V 1“C 

P + 


a 


A 

/'i, 




(2.49) 


6^97 
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¥ith. 3-^^^ determined from eqn. (2,49)? the emergent distri- 
bution for approximation is then, 


4(0, p) = I (%-~) 


T1=:0 


’ll 


1 


+ I (%— ) g„(-vi-^^) P^(p) 

n=o 


(1) 


2 ' ""n' o ' “n' 


(2.50) 


In P^ approximation, the two coefficients and a^ are 

obtained by the following conditions^ 

rl 




o 


( \ ap) 


. 3 ■ (-t-) «n(p) 

0=0 n=o 


5 

I 

n=o 


+ I (--2--) Sn a|i = 0 


■n 


and , 


•1 


j SjCiO { I ap) I 

0 ,i=o n=o 


2nPl ' 


+ I (— 2-b % (-'>¥') = 


(3) 


0 


n=o 


Simplification of these equations yields. 


s* 


P’ 3i 2('’P’) ^ Piptu) = 


- Pl2(-P5 


and 


^P^ P23(''P5 + ap' PpjtU) 


P 23 (-'-P 5 


f 5 ) (3) 

Solving these eqns, Si^ and a^ are 
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,( 3 ) 


‘^23 ^23 


rrr; 


and 


a 


(3) *^23^^ 


Pi2(-P^)-P..(vp)) 


^l2^^o 


1 


(3) . 


P23(-vp)) 




where v^ ^ is the largest positive root of g 4 (v) = Oj 
■with g^(^^) given hy eqn, (2.11) and p functions 
already defined by eqn, (2.37). The emergent angular 
flux for approximation is then. 


1 


’<'(0, p) = I a^.^^ t (-^-t^-) g (v.) P (p) 
3=0 n=o ^ 


+ I {-H-) h'f*’ <2.51) 


n=o 


Eqns, (2.48) and (2.51) are numerically compared with exact 
emergent distribution and the results are tabulated in 
Table 2,2. It is seen that the TPj^-approximation provides 
significant improvements over the conventional approxi- 
mation. However, the calculational results show that TP^ 
is not always better than TP^^ in the range of -1 h P' 1. 0. 

This is especially t rue for low values of c and p nearing 0. 
The exact solution )|/(0, p) is identically zero for p > 0, 

Any approximate solution will, however, produce oscillations 
of ^(0, p) for 0 < p <. 1 [44]. Of course, these oscillations 



Ta^ble 2^2as Emergent angular distribution for 
source-free Milne problem, c = 0.2 


1-i 

^1 


1-3 

TP3 

Exact 

-1.0 

1.2702 

0.1096x10^ 

3.4370 

0.1096x10^ 

0.1095x10 

-0.9 

1.1939 

0.9991 

2. 6665 

0.9943 

0.9936 

-0.8 

1.1177 

0.4998 

1.9743 

0.5037 

0.4957 

-0.7 

1.0415 

0.3332 

1.4003 

0.3416 

0.3295 

—0.6 

0.9653 

0.2499 

0.9343 

0.2591 

0.2462 

-0.5 

0.8891 

0.1999 

0.5662 

0.2068 

0.1961 

-0.4 

0.8129 

0.1666 

0.2857 

0.1688 

0.1625 

-0.3 

0.7367 

0.1428 

0.0827 

0.1385 

0.1384 

-0.2 

0.6605 

0.1249 

-0.0529 

0.1130 

0.1200 

-0.1 

0.5843 

0.1110 

-0.1315 

0.0909 

0.1053 

0.0 

0.5081 

0.0999 

-0.1631 

0.0718 

0.0924 

0.1 

0.4319 

0.0909 

-0.1580 

0.0555 


0,2 

0.3556 

0.0833 

-0.1263 

0.0421 


0.3 

0.2794 

0.0769 

-0.7820 

0.0319 


0.4 

0.2032 

0.0714 

-0.2389 

0.0254 


0.5 

0.1270 

0,0666 

0.0264 

0.0229 


.0.6 

0.0508 

0.0624 

' 0. 0626 

0.0249 


0.7 

0.0254 

0.0587 

0.0745 

0.0321 


0.8 

0.1016 

0.0554 

0.0518 

0.0449 


0.9 

0.1778 

0.0522 

-0.0155 

0.0638 


1.0 

0.2540 - 

-0.3753 

-0.1377 

-0.3353 
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Table 2.2b; j-imergent angular distribution fcr source-free 
Milne problem, c = 0,6. 


1-^ 

b 

TP^ 

P 

^5 

to 

Exact 

-1.0 

1.0552 

5.2195 

1.8602 

5.1451 

5.2065 

-0.9 

. 0.9918 

1.6168 

1.5596 

1.5915 

1.6056 

1 

O 

; ♦ 

CO 

0.9285 

1.0744 

1.2575 

1.0795 

1.0605 

-0.7 

0.8652 

0 . 8011 

1.0114 

0.8212 

0.7860 

-0.6 

0.8019 

0.6561 

0.7989 

0.6578 

0.6197 

-0.5 

0.7586 

0.5254 

0.6176 

0.5577 

0.5074 

-0.4 

0.6753 

0.4456 

0.4651 

0 . 4400 

0.4256 

-0.5 

0.6120 

0.5851 

0.5390 

0.5555 

, 0.5626 

-0.2 

0.5487 

0.5575 

0.2568 

0.2800 

0.5115 

-0.1 

0.4854 

0.2985 

0.1562 

0.2118 

0.2671 

O 

r : 

O 

0.4221 

0,2655 

0,0947 

0.1508 

0.2225 

0.1 

0,5587 

0.2571 

0.0499 

0.0977 


0.2 

0.2954 

0.2118 

0,0194 

0.0555 


0.5 

0.2521 

■ 0.1885 

0,0008 

0.0195 


0.4 

0.1688 

0.1660 

-0.0085 

-0.0055 


0.5 

0.1055 

0.1455 

-0.0104 

-0.0155 


0.6 

0.0422' 

0.1186 

-0.0078 

-0.0099 


0,7 

0.0211 

0.0890 

-0.0050 

0.0062 


0,8 

0.0844 

0.0481 

0.0016 

0.0512 


0,9 

0.1477 

-0.0227 

0.0056 

0.0487 


1.0 

0.2110 

-0.2146 

0.0007 

-0.0299 
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Table 2,2c% Emergent angular distribution for source- 




free Ki3.ne 

problem. 

c = 0.9. 


!-A 




LZZ.,.. 

TP3 

Exact 

-1.0 

0.6687 

0,8195 

0.7420 

0.7608 

0.8080 

-0.9 

0 . 6286 

0.7204 

0.6928 

0.6968 

0.7081 

-0.8 

0,5884 

0.6381 

0.6388 

0.6362 

0.6248 

-0.7 

0.5483 

0.5683 

0.5810 

0.5767 

0.5539 

—0 . 6 

0.5082 

0.5079 

0.5205 

0.5168 

0.4922 

-0.5 

0.4680 

0.4549 

0.4583 

0.4564 

0.4375 

-0.4 

0.4279 

0.4075 

0.3955 

0.3955 

0.3880 

-0.5 

0.3878 

0.3646 

0.3332 

0.3348 

0.3424 

-0.2 

0.3477 

0.3252 

0.2725 

0.2750 

0.2991 

-0.1 

0.3076 

0.2884 

0.2143 

0.2173 

0.2564 

O 

• 

o 

0.2674 

0.2537 

0.1598 

0.1627 

0.2073 

0,1 

0.2274 

0.2204 

0.1100 

0.1124 


0.2 

0.1872 

0.1879 

0.0660 

0.0677 


0.3 

0.1471 

0.1557 

0.0289 

0.0299 


0.4 

0.1070 

0.1233 

-0.0004 

0.0003 


0.5 

0.0668 

0.0902 

-0.0206 

-0.0200 


0,6 

0.0267 

0.0555 

-0.0308 

-0.0298 


0.7 

-0.0134 

0.0185 

-0.0299 

-0.0281 


0.8 

-0.0535 

-0.0218 

-0.0170 

-0.0142 


0.9 

-0.0936 

-0.0668 

0.0095 

0.0124 


H 

m 

o 

-0.1337 

-0.1186 - 

0.0500 

0.0516 
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are jin general ^ much less in approximation when com- 
pared with the Pj;j—approximation, Also, TP^ oscillation is 
found to be more pronounced than TP^ case. It can be thus 
concluded that oyer the whole range of p e (-1, 1),TP^ 
emergent distribution is better than TP^ result. Phis 
fact will also be reflected from leakage calculation to 
follow in the next section. 

2.5,2 leakage for Source-free Milne Problem; 


The leakage from the half space is defined as. 


L 



p p) dp 


(2.52) 


The numerical calculation of leakage, which is an integral 
parameter involving the emergent angular distribution, 
can also provide an overall comparison of • 'emergent angular 
distribution. In. conventional Pj^j- approximation and exact 
trea,tment the contribution of the integral in eqn. (- 2 , 52 ) 
for the range of p e (0, 1) is zero. This follows from 
I'larshak free-surface condition or exact boundary condition 
Because of the modification of Marshak conditions, this 
is not so in TPjj approximation. It is to be noted that the 
corresponding free-surfacc condition in TPj^ approximation 
is given by, 

rl 

(p + ajj) ,p(0, p) dp = 

■Jo 


0 



58 


Tfliich. implicitly in.corpora.tes higher order moments of 
angular density through Thus, the inward current 

from the vacuujn to the medium is not zero in this ca.S 3 
and hence the lea.kage in TP|^ approximation is given by 
e'qn,(2.52). The lea,ka,ge is given by, 


dO 


.-'O 


J -1 


\x Jp(0, p.) dp 


By use of eqn. (2,50) and (2.51) ^3 leakages are 


,( 1 ) 


= a 




(1) 




(2.53) 


= ap! S 3 _(vp 5 + ap^ - gi(vpb 


(3)^/^ „(3) 


(3---a^ ) + a 


(3) 


gl(v o' ■ J -^2. *^1^^1 


(vt) (2.54) 


For TP^ and TP^ approximations we get, 

+ Ap* gj_{Dp5 


(2.55) 

( 2 . 56 ) 


The exact leakage is [5] 


L = -.•2vJ X (-Vq)- (1-c) 


(2.57) 


Calculations based on eqns. (2,53)”(2.57) show that leakage 
improves in TPj^ approximations, in general. Table 2.3 
produces the results on leakage from the half-space. It is 



Talole 2.5? Leakage out of a source-free 


c 




TP, 

0.9 



0.2675 

0.2755 



0.2779 



0.2774 

0.8 

0.5405 

0.5611 

0.5721 

0.5711 

0.7 

0.5874 

0.4290 

0.4461 

0.4447 

0,6 

0.4221 

0.4915 

0.5126 

O..5II5 

0.5 

0.4495 

0.5552 

0.5754 

0.5764 

0.4 

0.4721 

0.6254 

0.6565 

0.6444 

0.5 

0.4914 

0.7066 

0.6964 

0.7206 

0.2 

O.5O8I 

0.8004 

0.7558 

0.8092 


half space. 


Exact 

0.2772 

0.5713 

0.4466 

0.5170 

0.5881 

0.6627 

0.7427 

0.8266 
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to bo noted that TP^ leakage is always better than TP^ 
leakage -values in spite of the fact that TP^ emergent 
distribution may not be bcbter than IP-] result for 
f. P- 1 0* This proves the fact TP^ oscillations of 
’l'(0, [j,) for \i > 0 is much less than TP^ oscillations. 

This is observed in numerical calculation of '^'(0, p) as 
shown in Table 2,2. 

2,5.3 Scalar Flux for c = 1: 

Earlier calcula.tions on source-free Milne problem 
have been done for c < 1, The special case c = 1 needs 
to be treated separately as demonstrated in the evaluation 
of a^jCc) fof c = 1. As a check on the usefulness of TP^^j 
approximation scheme in this case, the scalar flux <|)(x) 
calculated for P^ and TP^j approximations is compared with 
the exact va^lues. By integrating eqn. (2.38) with respect 
to l-t over (-1, 1) we get for TPjj approxims^tion, 

<4,(z) = ZqCI) -I- + I e ^ (2.58) 

In source-free Milne problem, the use of exact value of 

= 0,7104 ensures exact asymptotic flux for the case 

(IJ) 

of c = 1. The coefficients A^- ^ are evaluated from 
eqns. (2,25b) as usual. For conventional Pj^-approximation, 
the scalar flux is given by, 

^(x) = + X +' I e (2.59) 

® J ^ 
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(jj) _ , 

and a^. are detormiiiGd "by the usual Marshak "boundary 
conditions (2.24). Por example; the approximation 
is evaluated (there is no in this case) with 

the help of eqn. (2.24) fis. 


'1 


p ^ p + I dp = 


0 


= I = 0.6667 


By use of first two Marshak "boundary conditions^ z^^^ and 
a^ ^ in P^-approximation can "be determined, ^ is found 

to "be equal to 0,7051. Similar calculation yields the 
coefficients for P^ approxiiaation. The exact scalar flux 
is given "by 

rl 

(x) = X + z +j A(v) e“‘^'^'^dv (2,60) 

° v^_i 

I'lark [45] lias evaluated the transient integral numerically 
and the present exact results are taken from his calcula- 
tions. Ta"ble 2.4 compares the approximate solutions given 
by eqns, (2,58) and (2.59) with the exact results, TP^^ 
scale-r flux is superior to P^ flux for large values of x, 
I’P^ improves P^ calculations for all but very small values 
of X and is nearly as accurate as P^ results for all values 
of X consideredhere. This proves that a single tmnsient 
term in TP^ approximation gives fairly accurate transient 
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Talle 2.i 

1- ; Milne p: 

roblem scalar flux 

for c = 

1. 

X 


_ 

^3 

3 


Exact 

0.0 

0 , 6667 

0.7104 

0.6023 

0.6072 

0.6108 

0.5773 

0.2 

0.8667 

0.9104 

0 . 8360 

0.8408 

0.8502 

0.8495 

0.4 

1.0667 

1.1104 

1.0585 

1.0635 

1.0725 

1.0731 

0.6 

1.2667 

1.3104 

1.2737 

1.2788 

1.2853 

1.2858 

o 

• 

CD 

1.4667 

1.5104 

1.4839 

1.4891 

1.4930 

1.4935 

1.0 

1.6667 

1.7104 

1.6908 

1.6960 

1.6977 

1.6985 

1.5 

2.1667 

2.2104 

2.1998 

2.2050 

2.2035 

2.2051 

2.0 

2.6667 

2.7104 

2.7031 

2.7084 

2.7058 

2.7079 
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flux representation in this case. It is observed, in 
general, that is superior to its Pj^ counterpart for 
large x* The reason for this in obvious 5 while TPjj 
approximation, through a- formulation, constructs an exont 
asjTuptotic flux in this case, the P^^- approximation is a 
total theory (not as3rmptotic and transient separately). 

2.5.4 Constant Source Milne Problem; 


As has been pointed out earlier, this problem 
explicitly uses the value of in TP^^ calculation. The 
exact angular flux in this case is [5] 


q -x/v 

4;(x, la) = TTPtn- + a^ 4 - ® ° *^0+^^^ 


7_ 1- a^ , y 

( 1 -cj 0+ 

rl 

o 


+ A(v) ‘f'v(u) 


( 2 . 61 ) 


where, 


a 


0+ 




( 1 -c) CVq X (vq) 


A(v) = + v) X (-v) g(c,v) 


g(c,v) = [( 1 -cvtanh” v) + 


-1 o 2 2 2 T 

1 %2 , -n: c V -,-1 


and is a constant isotropic source. 

The exact leakage is computed from 
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j \x i)(0, ij.) &\x 


wh.Gi'e ip(0, ij.) is obtained from equation (2.61). This gives, 
2q 

L = ("o “ 


rl 2 

1(c) = I \ [1 + ] qg(c, \i) d]X 

Jo (1-u ) 


where, 

rl 

[1 + 

'O (1-.^") 

Values of 1(c) are tabulated by Case, de Hoffmann and 
Placzek [23]. The angular flux in approximation is 

1. 


il) (x, p) = + I e ^ 




z 


1 (%— ) 


21=0 


In this case, coefficients a(^^ are determined by usual 
I'larshak boundary conditions (2,24), Thus for P2_ approxi- 
mation we get, 


or 


i.e. 


-^0 n =0 


0 


= ( 




For Pj approximation one gets. 
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a 


(3) 


2TI:::ht 




Pia'^i) P23("P^) - Pipt-P^ p,^(v,) 


and 


12^ o ^ ^^23^''!' 

( 2 . 62 ) 


a(3) 

1 




(3)^ o l .( 3 ) 


^12^^!^ ^23 ^ " Pl2^^o )P23^^1^ 

(2,63) 

wlioro Pjj^^ are defined by equation (2.37). 

The TPjjj flux is given as 


* (x, n) = 7T~rv- + 


T1~S)- %+ 

+ i ' aW ? ,.^+1 

3 ^ 


i_ gyv.)py,) 


n=o 


In this case, it is to be noted that the as 3 nnptotic coeffi- 

fn) 

cient is not equal to the exact asymptotic coefficient 

(equality holds for only source-free Milne problem on which 
the formulation of aj^ is based). Applying the modified free- 
surfacc boundary condition for TP^, one gets 


r 


1 


1 


(p -I- a;^) C ^rsi' ^oi^ '^0+^^'^^ = 


or 


*^0 1 , "''O 


(1-g 1 2 1 Tl^c) 0+ 

rl Cv 


(1) _h_q 


+ 

0+ 1 


0 


0 2(V^-.1J.) 


a 2(v^-p) 


dp = 0 


0 


dp 
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1 . e . , n. C 1 ) = - ] 

'0 4- CV ^ — 

- 2 ° [(''o+“i) In - 1] 

Similarly for TP^ approximation and Ap^ 


(2.64) 


0 ^( 3 ) _ ^^0 3 ^ ^1 2 1 ^ *^^ 3 % l ^ ^ 23 ^'^ l ^ ^^ 3 *^ 

° 2 ^12^^ l^'^'^3%l(^l)] ^31"P23^’^1^ ^10 

(2.65) 


Ap) = „ 

"2 I^Pi2^^1^ ■*■ “3p3i(^l)] ^3i“p23^'’l^ ^10 

(2.66) 

whore , 


I 


10 


I 


31 


+ ctj 

51 (3) _ 3i-(-l) 


a-, and a-r in tlio ahove ecuations are obtained from 
1 

eqns. (2,29) and (2,35)» 1^5 calculation, a^(c) is 

obtained from eqn. (2.44). It is observed from numerical 
calculation that the asymptotic coefficient in TP^jj. calcu- 
lation is superior to the Pj^ values, in general. Table 2.5 
compares tho values of asymptotic coefficients in TP 
approximation with exact and P^^ coefficients when I = 5. 
Vath the knowledge of coefficients ap^ and A^^ TPj^ 
loakago can now be easily evaluated as follows. 


I: 



i 
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Table 2,5 t The asymptotic coefficient for N = 5 
for a h-alf-space of constant source 
strengtli. 


c 

a(5) 

P 

0 + 

O' . 

04- 

0*9 

-5.6485 

-3.6443 

-3.6447 

0.8 

-2.0780 

-2.0723 

-2.0758 

0.7 

-1.3677 

1.3483 

-1.3550 

0.6 

-0,9554 

-0.3929 

-0.8941 

0.5 

-0,6494 

-0.5773 

-0.5505 

0,4 

-0.4582 

-0.3653 ^ 

-0.2806 

0.5 

-0.5309 

-0.2279 

-0.0896 

0.2 

-0.2458 

-0.1480 

-0.0073 
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Iioakagc 


(Tl) 


r-L 

j-1 


|i) d!,i 


r V 


( 1 ) 


.1 '-rr--5T + an 

si{''o) 


PPjj Leakage 


(TN) 


a 


(10 

0 + 


®l(»„) + 




2 


. x '- 7 ’> 


In conventional approximation, these expressions change to 

a.nd „ 




L 


W = I' af)g^(.W) 

i-1 ? J -i- J 


Loako.go V 3 . 1 UGS with a constant source q = y"5(l-c)' for exact 
solution, TPjvj and conventional Pjj- approximations are 
produced in Tahlo 2.6, Numerical results show that TPjj- 
approximation does not necessarily improve leakage for high 
absorption (i.o., small c) cases. The reason may be 
attributed to poor transient flux representation for high 
absorptions in TPj^ approximations. Thus although the 
asymptotic flux is improved in TPjj approximation (Table 
2,5 shows this), the advantage is suppressed by a corresponding 
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Table 2.6; 


leakag© Trom a half— space of consfani} 


source strength. 


c 



TP^ 


TP3 


TP 3 

Exact 

0.9 

0.7525 

0.6968 

0.7179 

0.7162 

0.7160 

0.7153 

0.7147 

0.3 

0.6595 

0.5S52 

0.6416 

0 . 6385 

0.6391 

0.6371 

0.6372 

0.7 

0.6126 

0.5085 

0.5930 

0.5862 

0.5900 

0.5843 

0 . 5878 

0.6 

0,5779 

0.4254 

0.5573 

0.5467 

0.5540 

0.5400 

0.5513 

0.5 

0.5505 

0.5576 

0.5291 

0.5250 

0.5255 

0.4985 

0.5227 

0.4 

0.5279 

0.2447 

0.5060 

0.5518 

0.5021 

0.1-612 

0.4990 

0.5 

0.5086 

0.1593 

0.4864 

0.6950 

0.4823 

0.4313 

0.4791 

0.2 

0.4918 

0.1154 

0.4694 

0.^323 

0.465 2 

0.4097 

0.4616 



70 


worsening of the transient flux for low values of c. Since 
loo.kagc is an integral parameter involving angular flux, 
bhis fact is reflected in the leakage results. The failure 
of TPj^-approximation for small values of c in this case 
calls for a hotter transient flux representation in the 
approximation scheme. However, it is felt that this may 
not he realized in the ; .framework of approximation 
and a different representation for the transients may he 
necessary to achieve any tangible result. The other 
possible remedy may he to separate two class of transport 
problems [46j - one more close to source-free Milne 
problem and the other reasonably approximate to constant 
source Milne problem. In that case, TPj^ approximation 
needs to bo formulated separately for constant source 
typo problems. A detailed discussion along these lines 
has been given j.n Chapter IV, 


2,5.5 Critical Problem: 


The present sot of calculations involving a free- 
surfaco culminates in the study of critical problem in 
slab goomotry. The exact free surface condition ‘ for a 
slab of half-thickness b is [5] 


b-A^ , 


.b„A 


fl .. ♦ (a) + a-'' = ° 

+ A(v)[e ° 

TO 




(2.67) 
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iiqii. (2.67) is the starting equation from which the 
cri'cical half-thickness o.ncl corresponding flux distrihu- 
tion are determined in exact analysis. In conventional 
Pjj approximation the free— surface condition is 


(h) 

c'’o r ,2n-:-l>, , (u), , , 


b„Ay’' H 


n=o 


‘n 


+ 


-b /v 

f\ f 


(I) 


IT 




4- 


Z'aW [ 


n=o 

I 


^ l (-2-) g^(v ) P (^) 
n=o ^ 


-b„/v. I 




0 


n=o 


}i >_ 0 (2-68“-J) 


The corresponding TP^x free-surface condition is 


■N 

b_/v “h /v 


1 


(If) b /v. y (^^1) g (v .)p (u) 
+ r A, [e ^ ^ nio 2 3^ n^l"^ 


'■ ,n ,2n+l. 


U -1 

3 =roots 

+ e ^ I (-I)"" (^“) g^(v JP^(p)] . = 0 

n=o * 

p 0 (2,68^ 

Because of the homogeneous nature of eqns.(2,67) and (2,68) 
assnnptotic coefficient has been put equal to unity by 
normalization. In eqn* ( 2 , 6 %^ the unknowns {b, 
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are d etermined by application of conventional Marshak 

boundary conditions while the modified free-surface 

conditioiis as given in egn. (2.25) are used to evaluate 

the unknown coefiicients of eqn, (2.68'). For approximation, 

one gets for the critical thickness b^^^ 

c 


or. 


4- G 


I (-1)“ (5Sii) gjj(k^’)P„(ii)] du = 0 

n=o ii- o ij. 

( 1 ) 


b /v 
c' 0 


( 1 ) 


[ ^ g^(v )] + e ^ 






= 0 


By use of = ik^ and noting that 


b 

. c' 0 


b b 

= cos ^ - i sin 

i'ju _ is- ... 

0 0 


\I e get , 


or 

where 

length, 


cos 


,( 1 ) 




0 


+ k^^^(l-c) sin — = 0 

^0 

.( 1 ) 




k 


( 1 ) 

D 

,( 1 ) 


l/lf^Tc-lT and linear extrapolation 
2 / 3 . 


Similarly, the critical half-thickness for TP^ approximation 
is obtained from the equation 
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(|i + a^) [e ° <I>q+(!J.) + 


-b /v 

® "" ° '^o»-OjO] <3p. = 0 


B 3 ^ use of = ik^ , the expression 

^./v. 


n' ~n -b / V 

® H- e <= » <,^_(^) 


(2.70) 


simplifies to 


ck2 
0 


b ck u 

2 oos i.- - 

0 k„ + |i 


3in 


k 


fheii the condition (2,70) is written as 


i 


■■'1 ck^ h ck p "b 

(p-i-a, ) (• — « — fj— cos ^ 'A —- -'5 sin ~) £jp_ 

0 k^+p^- 1^0 kj + p2 k^ 


= 0 


By use of from eqn. (2.30) it follows, 

p 

b ]r 
‘^0 ir- = ‘T*-' 


(Tl) 


1. e. 


k^ tan“^ (k^/A) 


(2.71) 


by their exact lvalues 


which is the end-point result [5]. E(in. (2,71) clearly 
shows that replacement of different Pj^-approximation 
parameters like 
(v^, X f z^) is indeed possible by a modification of free- 

surface conditions. 

In TP^-approximation, the two equations in the 


(3) 

unknowns b^ and are 
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A 

j ( [J'+OC^ ) 






° VW -e . (,) 


f^)r r /2n+l, , , 

I_ (—2-) gjj(vj^) P^(n) 


+ "^1 L ^ 


n=o 


^ ® X Sn'^) = 0 

( 2 . 72 ) 


and, 


fi 

0 


b / V 


° ° + e 


C 0 


d^i 


(3) 

-b^/vi 3 


lo j 

+ Ai- [C <= “" Jo (- 2 - 1 ) ,, 


-f e 


ti-o (-1)^ i—'^) P^^(lJ.)P^(^)dii] 

(2.73) 


= 0 


I'/e already have, 


" 13 . /v^ 
<!'o+ (l-i) -I- e ° 


ck.^ b_ ch u 

rr"""^' ir ■ :t" 

ko + li o k^+|i 


2 ^ 2 ' k 


The integrals o:f the type 






n 


-'o (k" + p^) 


are evaluated by ^ke following recurrence relation [43] 
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J(n) = 


' 1 v2 


0 (icr + \x^) 


2 ; = iT-i - < J(n-2) 


irith, J(0) = tan"-^ (^) and J(l) = i 

■"0 ^-^0 ^ 

Thus, eqn. (2.72) and (2,75) simplify to 

b 
i 

^0 

b 

- sin [ck J(2) + ck a^J(l)] 

0 ^ J 

+ [e ^ Ptp (vT ) + e ^ 


In 


2 

k^+1 

o 

* 'T- 


cos Tr [ck^ J( 1 ) .|. ck^ J(0)] 

n ^ J O 


“H OC^ 0 




p-| j(v^) + Oij ® Pp^(~vp)] = 0 

(2.74) 


and , 


ck^ COG ^ [ I J(5) - I J(l)] ■ 

0 

b ^ 

- ck sin ]“' (^- J(4) - ^ J(2) ) 

O 

/V h'/v^ "’b^/^1 

+ Ap> [e = ^ P 23 <'’!>+ ® ^ P 23 (-''l)] = 0 

(2.75) 

where p are given by eqn. (2.37). Bqn. (2.74) and (2.75) 

" (3) 

are two homogeneous eqns. involving b^ and . Eliminating 
j^(3) these equations it follows, 
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{GOB ck^ J(1) + ck^ ct^J(O) - sin ck J(2) 

0 ^ 


oko sin ajJ(l)}{ P 2 j(vj^) + e 


-■b /, 


^23 ^ 


- { ck|^ cos J(3) - |. J(l)) -- ck^ sin j“(| J(4)“ |j(2))l 


^ ^5 

O ^ ; ~ y d \,X J J - C 

0 ~ ^ 0 

^c/\ -6^/v- 


K/ 


+ [e ° ^ + e ^ = ° (2.76) 


-6. A 


Tile solution of "this ■transcendental ei^ua'tion gi'ves "tlie 
critical half thickness ' h^ for TP^ approximation. For 
the solution of (2.76), Newton-Eaphson method in employed. 
a^(c) in obtained from eq.n. (2.36). The critical half - 
thickness for P^ approximation is obtained by solving a 
similar tranGcondental equation. Table 2.7 gives the 
critical hr.lf-thicknessos which sho'w that the TP^ approxi- 
ma.tion imprcu-es the end-point result (i.e. TP 2 _ result) 
further. 

2 . 5.6 Interface Problem; 


For an analysis of interface problem by TP jj method, 
the two half-space Milne problem is considered as the 
standard reference problem. The angular flux i|)(x, p.) of 
the two half-space Milne problem satisfies the following 

equation. 



77 


Talile 

2,7; Critical h.al: 

mean, free pa' 

f-thiclaiess b 

c 

■'ull. 

in units 

of 

c 

^1 


^3 

TP^ 

Exact 

1.02 

5*7519 

5.6655 

5.6711 

5.6655 

5.6655 

1.05 

3.4054 

3.3002 

3.3066 

3.3002 

3.3002 

1.10 

2.2207 

2.1134 

2.1213 

2.1134 

2.1134 

1,20 

1.4125 

1.2893 

1.3020 

1.2895 

1.2893 

1,40 

0.8531 

0.7384 

0.7577 ■ 

0.7364 

0.7366 

1,60 

0.6269 

0.5152 

0.5384 

0.5113 

0.5120 

1.80 

0.4965 

0.3932 

0.4182 

0.3882 

0.3887 

2.00 

0.4121 

0.3163 

0.3421 

0.3110 

0.3108 
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l-i II' (x, \x) + la) = ° 


2 j ^ ’^'(x, ji' 


) du' X 0 


I H 

J -1 


X, ii' 


■i*rii:]i the following boundary conditions 


lim \p(x, |i) 

X ■+■“ 

lim 4;(x, |.i) 

X 4-- -00 




dji’ x > 0 


(2.77) 


and the interface condition 


;(i(o+, lO = ^(o--, ii) 


where c and c are number of secondaries per primary for 
the loft hand and right hand media respectively. 

The exact solution of eqn, (2,77) is [5] 

x/v fO -x/v 


I|j(x, IJ.) 


.a„ e 
0- 


.0 ~x/v 

X < 0 (2.78a) 


'iix, ti) = e ° >'o„(^) + %+ ® 
A -x/v_ „ 




+ ( A(v) e , x>0 (2*78b) 


The coefficients a^_,, and A(v) are given as [5] 
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where , ) 

and> 

X and g' 
defined 


and, 


a. 


x(Vo) (v + V ) 


a. 


x(“Vq) 


2 c 


X(-V^) cv^ (v^ + v^) 


A ( V ) 


, V > 0 


A ( V ) 


v‘c(l-c) (v^- x(-Vq) X (»v) 


(Vq + v) I(v ) X(-v) 


vc(l~-c)v^ (v^- v^) x(-Vq) X(v)(v^-v) 
(Vq - v^-) ll(v) X(v) 


V < 0 


function is defined as 

x(k '0 = x(z) x(-z) 


N(-) = "7—r 

g(c, ) 

-functions in these equations have already been 
, In the TPjj formalism, the solution (2.78) becomes. 


li-Cx, |i) 


i (x, [i) 


x/v 


-aL« ♦oJp) - 


IT 


I Sn (-^i) X < 0 (2.79a) 

n=o 


x/v 

® ' <!>o- 


° ” (i*) + 

K 




I'l'W I (^-b gn(vj)P„(|i). 


n=o 


X >0 (2.79b) 
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Where prime in the j-summation indicates that only posi- 
tive roots are taken, 'the interface condition then is 

- i'aW ] (2a+i) g^(. ) J { ) 

0 n=o j n 

2 n=o 

(2.S0) 

Multiplying both sides of eqn.(2.80) hy P^(p) and inte- 
grating with respect to p- between -1 and 1, it follows 


-a 


(ir) 


s„(-vj - I aW 


3 


= ^(-V„) g^(vp +lh<®)g^(vp 


3 


n = 0,1,2, ..N (2.81) 

Eqn. (2,81) is just of the same form as the interface 
condition used in conventional approximation. In 
approximation, we wish to have an exact transport theory 
asymptotic flux for the Milne problem treating the 
adjoining half-space as purely absorbing (c =0), Thus a 
parameter 6i.j(c) is introduced in the first equation of the 
set of eqns, (2.81) as 


-a 


(K) 




(N) 


(tT) 




1 *U < %+ 




(2.82) 


Note that 6jj(c) is introduced as a multiplying factor 
in til© non-source part of the solution. The eqn. (2,82) 
alongwith the following eqns 


(l^) 


„(N) , v\(l) I . 

*a„ . ^ ^ gv,(-vJ 


’0"-' 'o 




J "n' ''j 


Sn(-Vo) + 


'’■o7 En("o^ I At®) ant''j) 


n = 1,2,3, ...N (2,83) 

constitute the interface conditions for evaluation of the 
coefficients. 

Determination of 6^(c) for TP-, and TP^ Approsimations ; 




6jj(c) is determined "by the condition that exact 
asymptotic flux he obtained for the right side medium 
treating the adjoining lialf~space as a purely absorbing 
one i*Q. c = 0. In that case, takes its exact value 


i. e. , 


a 


(n) 


a. 


0 + 


■2z’ /v' 
0 0 


Ot “ ‘ 0+ 

In approximation, it follows from eqn. (2.82) 

(1) 


■a. 


l“f" C-i Q* 


1 0 + 


(2.84) 


when c s: 0, the asjmiptotic flux in the left hand side 
medium is 7<0ro. Therefore, we get, 


a 


( 1 ) 


0 


The condition (2.84) then becomes 


1 + 6n a 


1 “0+ 


i* 01 


2z /v 
o' 0 


(2.85) 


0 + 


0 
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For IPj-approximation, from eqns. (2.82) and (2.83) the 
interface conditions are as follows.. 

_a(3) .(3) ■ ,-(3) -(3), 

0- H ■ = 1 + 6^ { > 

- a: (-V ) - A-f^^ (-\) ) 

0- ®n^ O'' •'^l 1'' 


gnC-Vo) + 


1 ^n' ''1' 


n = 1,2,5 


when c = 0, ap^ = 0 and a'lP 

0“r 


’^■0+ ® 


-2z /Z 

O 0 


Then, 6^ is determined from the following equations 

-Ap^ = 1 + 6. { + Ap) } (2.86) 

-Ap^ g^(-vi)|^^^ = + V g^(“o) 

+ n = 1.2 (2.87) 

In eqn, (2.87), 1° 1® evaluated at c = 0 and 

it is denoted as Sn("^l)lc=o 

From 0qns.(2,86) and (2.87), it follows 

Bjj2 ®22 

1 + — g-- J— 


“11 12 

®21 ®22 _ 

®21 ®E2 -2^0'^c 

- e 

ni ^12 ^ 


( 2 . 88 ) 


®21 ^22 
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where , 


and, 


in 

■'■11 

O 

II 

a 

> 

II 

’ *21 - 

•“g2^^l) 

*12 

If 

f 

, E 22 = 

-g2(\) 





*R1 

= (-e ° ° 

“1) SiC^o^ 



-2z. / V 


%2 

= (-e ° ° 

+ *) g2^~o) 



c=o 


The values of 6-^ and 6^ as obtained from eqns. (2,85) and 
(2,88) will be used for any problem involving interface. 

Puel-Iioderator Cell Problem: 


An infinite homogeneous array of 2w cm wide 
fuel plates separated by a moderator with a centre-line 
spacing of W cm is considered. The disadvantage factor 
'd' in such an array is defined as the ratio of average 
flux in the moderator to the average flux in the fuel i.e,, 


d 



C(x) dx 


'l>(x) dx 


(2,89) 


where <j! (x) is the scalar flux for the fuel and (x) is that 
for the moderator. All the quantities under the bar are 


used in case of moderator. The disadvantage factor Vd' 
will he calculated by TPj^ formalism and the results will 
be compared with values obtained by conventional P^ . 
approximation and exact solution. 
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In thr^ TP^ approximtion, noting that angular 
fliucen both in the fuel and moderator regions are symmetric 
abou'G 'clioi.}; ;iiid~pla,ne 3 the angular fluzes are given as, 
the origin at the centre of the fuel) 


li) 


\x) 


C-r 


{ 

0+ ^ 


cv. 


Zz/' 


CV 


2(%+li) 




2(Vl^) 


-Zx/\ 


( 2 , 90 ) 


c V, 


'z(z-W)/v^ 


2(Vo+lJ.) 




c V. 


Z (z-¥ ) /v 


2(v^-|l) 


° ,.--S 


(2.91) 


(1-0) 

A spatially constant source has been taken in the moderator 
region with no source in the fuel. The factors z and f 
arc incorporated in eqns. (2.90) and (2.91) because distances 
in this case arc measured in cms. The interface conditions 
at the interface z .= w are obtained from eq_ns. (2,90) 
and (2,91) as 


and, 


a <9 ooehS 

0+ ^0 X=¥ 


( 1 ) 


6, a^^^ cosh 
1 0+ 


g^(vQ) sinh 


zx I 


“0+ '=>1^ 'O' 


V fx=W (1-c) 
0 

(2,92) 


g ("v ) sinh 

®-o+ ®1^ oh “ 


(2.95) 


z=w 
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It is to be noted that with the introduction of the 

parameter the moderator scalar flux is given by the 

right hand side of eqn. (2.92). Thus the scalar flux at 

the interface is concinuous -» this is at variance with 

Pomraning s aS3niiptotic L/heory [26] or the usual as 3 niiptotic 

diffusion theory where s. discontinuity is allowed to 

improve the overall result.. In addition, the present 

formulation by its construction yields a better angular flux 
fl’) — d') 

a^:j'/and are given from eqns. (2.92) and (2.95) as 



(1-c) sinh cosh ~ 


+ g^(v„) sinh — cosh > (2.94) 

JL J- 0 


a 


(1) _ 


g- (v ) sinh “ 

o^ Vq 


0 + 


( 2 . 


(1-c) cosh ™ 


+ 5-, gn (v ) sinh “ cosh ^ > 

1 i. 0 V v' 

0 


The disadvantage factor is obtained from eqn. (2.89) as 


i 

(l-c) 


+ 6 ^( 1 ) ._lo_ sinh 

^ V ^-o*- (W-w) z \ 


w£ 


Ji„. a sinh ~~7 

w 0+ 2; V 


95) 
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By uso of and from oqns. (2.94) and (2.95), this 

expreoDion aimplifies to 


a = VS ooth 5 '. + 6 , .-.i- ilr 4 ^_ 

"o -0 1 (iJTj 


X { (¥-w) Doth ~sv-vl _1 ) 


(2.96) 


Til© convontlonal “approximation (3isa(3vantage factor 
follow© from eqn. (2,56) ly replacing v^, by their 
value© and by putting i't follows for 

Pj^-approximation, 


a = ooth 


EW ^ W (l-c)E 

(1-c) r 


0 


X { (l'/-w) coth - 1 > (2,97) 


In TP^-iapproximo.tion; the angular flux in the fuel and 
moderator sides are given as 

(x, li) = 


,_s CV^ r. -Ex/v 

a ( 5 ){ P ^ 0 ^ 3 


0+ 2i\n) 


+ aP\' i (“f^-) gj^(-v^) P^(p) e 


aTv^HT 


£X 




'1 


n=o 


+ I (i|.t±.) gj^(v^) P^([i) e 
n=o 


- Ex/ > 


} (2 


.98) 
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'P(x, |.l) 


(3) r "o 




z(x-¥)/7 


2(vQ+!a) 


c V 

-1- 

2(VQ~ia) 


“-'e( x~W ) /v 


A 


(5) ( ? /2n+lY / - 

- ^ 1 (-p-‘-) g^(-v,)P (p.) e 


1 ' V 2 ' 

n=o 


+ n-^r--) «n(vi)Pn(l^) ® 

n=o ’ 


-I (x-W ) /v . 
e ^}+ 

(1-c) 


( 2 . 99 ) 

Tho intorfocG conditions for TP^-approximation are olDtained 
at X = w 8.0 V 


a 


GO all -^r- + cosh r~ 


0 + 


' '■*3 

0 


6- < K(3) cosh + If' cosh 1 

, 0. 1 

+ — -? — (2.100) 


(1-c) 


a 


(3) 

0 + 


gj (v^) sinh + Ap^ SlA^) Slnh “ 


Ew 




+ Ap^ g^(v^) sinh -^£— > 


( 2 . 101 ) 
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(3) 


'o) «;■ Sjf"!) =osll § 


Ew 

1 


'^0+^ S2^^o' 


i:,(¥~w) 

‘v 


+ Ap) gj(~^) cosh 


'1 


( 2 . 102 ) 


'^oV ®3('’o) ■''“'1 i;- + 4^^ s^(\) slah ^- 

= - { a^:j'.^ SjG’g) sinh Sfcjll. 

^0 

+ Ap^ C''.(v. ) sinh lilTll } (2.105) 


Usin^i' 65(c) from oqn. (2.88), eqns. (2.100)-(2.105) can be 
solved for tho coefficients aP^, aP^, Ap^ and iPP The 
TP^ dif^advantage factor is then calculated from eqn.(2,89). 
Tho conventional result follows by replacing the asymptotic 
roots by their P5 values and by putting 6^ = 1. In order to 
compare the TPjj results with conventional results, tr^o 
fuel-moderator cells given by Murray [47] and Ifeinberg 
and V/igner [44] have been considered. The results are shown 
in Tables 2.8 and 2.9 respectively. It is found that 
TP|j approximation improves the result significantly. 
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XciTdIo 2* St C©!! fcic'toi' foir ppot)l 6 iii 1 [ 4 .Y J 


w 

= 0.635 cm 

w 

= 15.240 

cm 

B 

= 0.3970 

cm 

‘i 

a 

= 0.398 

-1 

cm 

t 

= 0.3748 

-1 

cm 


= 0.398 

cm“^ 

a 

= 3.31 X 

10 cm 

i' 

a 

CM 

K^ 

• 

0 

II 

cm 


Disadvantage 

factor 

Pi 

2.213 

pj theory 

2.405 

Exact 

.. 

TP^ theory 

2.330 

theory 

2.458 
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'J'nhlQ 2,9:; Cell disadvantage factor for 
Problem 2 [ 44 ]. 


¥ 


0.9294 

cm 

¥ 

= 29.3104 cm 



s: 

0.3951 

cm 


= 0.404819 cm' 

-1 

s 


0.3951 

cm 

^t 

= 0.404819 cm" 

-1 

z 

a 

=5 

0.3796 

-1 

cm 

I 

a 

= 3.81 X lO"^'' 

-1 

cm 


Disadvantage Percentage 



factor 

error 

theory 

5 .. 201 

5.30 

theory 

5.462 

0.51 

Exact 

5.490 

- 

TPi theory 

5.420 

1.27 

TP 2 "thoory 

5.482 

0.14 



Cli/IPTER III 


ANALYSIS OP PUTTIE gi^B PROBLEM 

Th© work of the present chapter concerns the 
application of the singular eigenfunction method for 
monoenergetic neutron transport in a finite slab. Prom 
what has been outlined in Chapter I, it is seen that the 
problem of neutron transport in a finite slab reduces to 
the solution of two coupled Predholm integral equations. 
Previous works in this area were mainly concerned with 
numerical solution of these integral equations. The 
present analysis also starts from the integral equations 
but differs in the important respect that the two Integral 
equations involving the discrete and continuiun coefficients 
are combined together to form a single Predholm integral 
equation. Function anal^'^tic study of this integial eqn. 
is the central theme of the present chapter. Per the finite 
Slab, the critical' problem is posed as an inverse eigenvalue 
problem where the eigenvalue is given and a parametric 
variation of the operator has to be detenained. When the 
eigenvalue is unity^ it defines the crltieal, slab psfoblesi 
and the corresponding inverse' parameter is the 
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half-tnickness. it is proved th^+ 

^ solution of finite 

slab critical problem is unioue a ^ 

unique, a good estimate of the 

inverse parameter (i p tho ^ • 4 .. . 

' • • “i-itxoal thloknese) is also 

obtained by applloatlon of bounds of s^ , 

ounas of spectral radius of 

Itnear poaitivs operator. This involves the solution 

of a transcendental equation which Is achieved by the 

method Of bisection. At the end of the chapter, a short 

discussion on uniqueness of solution of slab albedo problem 
is also given. 


-Eq uati ons o f Case/s Theory: 

We consider neutron transport in a slab of thick- 
ness 2b whloh is centered at the origin. In the case of 
slab albedo problem, a parallel beam of neutron is incident 
on the surface x = -b, while the absence of incident 
neutrons defines the critical problem. Mathematical repre- 
sent.'.'clon of the problem is given by the following equations 
When the scattering is isotropic 

1-^) = ^ ^( 2 :, p’) dp’ (3.1) 

'f’(~b, p) a ^(p-p^), p, p^ > 0 for albedo problem 
w 0> p >0 for critical problem 

?Ji(b, p) a O p < 0 

(% 2 ) 
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In Chapter II, the general solution of eqn. (3.1) 
in alroady given and discussed. In the following, we 
adopt the methods of McCormick and Mendelsoii [17] to 


obtain the coupled Fredholm equations referred to earlier. 
Thus the solution of (3.1) is written as, 

-x/v^ x/v 

0 - 


p) 


a e 
0 + 




x/v 


+ \ A(v) e <}>y (p) dv 


-1 


where the asymptotic roots i (-1, 1) are the solution 
of the eqn. 


Cv - v+l 

T ^ ‘v^r 


and the transient roots v e (-1, 1). The eigenmodes 
**’v, given by 






“r ^ vTu ^ 6(v-i^) 


T v-p 

where P denotes the principal value and A(v) is defined as 


^(v) 


c ^ taiih"" ^ 


To allow integration in the transient term between 0 to 1 
only, the angular density <J>(x, p) is written in the form 


° + - a - ^ V <^> 


-x/v 


dv 


'i>iXy p) « ® 

.1 

+ , A(V) ty(P') ® 

Jo 

A(-v). <^^v{p) ^ 

" J Q-- , :hhi:hK 


(3.3) 
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By applio' tion of bounOorv condition (5.2) m equation (5.3), 

bli© i oHowiii;',', SQ lift b ions ©.rs obtained for the coefficients 
aQ.|, anc'l ( v ) , 

6(p.-P-q) = 8 ” ^o- ® ° 

+ ( A(v ) "Jy (^i) d V 

Jo 

+ [ A(-v) (j) (ii) dv > 0 (3*4) 

-b/v b/v 

0 *= e e 

rl -b/v 

+ , A(v) ({) (li) e dv 

■J Q 

+ A(-v) «^(ij.) dv, u < 0 (3.5) 

-'O 

Changing the variable p. to -p in eqn. (3.5), gives, 

-b/v b/v 

® = %+ ® ° + ^0- ® 

+ f A(v) If (p) dv 

^0 

+ ' A(~v) i})^(p) dv, p > 0 (3.6) 

Jo 

Adding and subtracting eqns. (3.4) and (3.6) it follows, 

\ 

6(p~Pq) a (()q^(p) 

. . -Jb/yo i V : 

+ ® ° ^^•o- - ^"0+^ 
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^*1 

(A(v) _+ A(-v)) av 

ri , 

■h (A(™v) + A(v)) dv (3.7) 


Eqn. (3.7) is rewritten 

^ \ ~ 213 / V 

o(ii-ii^) = ♦^^(|,,) _+ b^ e ° 

(1 

+ Sj.( '') %(]i) dv 

A - 2b/v 

± . 53 .(v) e (]i) dv , la > 0 (5.8) 

uO - “V 


where the new coefficients b, and B,(v) are defined as 

i "t 

b/v 

5^4. * (ao4. ± a-Q^) e ° (3.9) 

B^(v) = (A(v)+A(~v)) (3.10) 


The coefficients b_j_ and B_|_(v) of eqn. (3.8) are evaluated 
by lialf-range orthogonality relations. It is known that 
the functions <1 q^.(P') o-nd 4)^ (p.) , 0 <_ v < 1 form a complete 

set for functions defined in the range 0 p, < 1. Also the 

functions bQ.j,(!-i') and ('>) are orthogonal with respect to 
a weight function W(p) = . yCt^) where. 


y(i^) 



2 (l*-c) 


tr 

Pi. 


(-p.) ( 




Thus, eqn. (3.8) is multiplied by V/(p). ^q^Cp-) and ¥(p) . 4 *^ (p) 
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and then integrated with respect to p between 0 and 1 to 
obts.in the following eouations for evaluation of b and 

ri 

6(i^-iio) w(|.) a;, ^ 


2Vv, 


+ e 


'+ 


1 


^(l^) dp. 


rl rl 

+\ B^(v ) , ¥(p) ^ (p) (p) dp dv 

■-'o • ' -j 0 


rl 

±'i B (v) e 

Jo -■^- 


■^0+ 
•2l?/v„ rl 


W(p) $Q^(p) t!>_y(M') dp d 

(3.11) 

fl rl 

\ 6(p-Pq) W(p) ^^(p) dp = b^i ¥(p) ^^(p) ,j^_^(p) dp 

JO ) A 


■2b/v, 


i ^ 

1 
0 


+ 


K (h) '!'o_(P-) dp 


+ ( B (v) ( w(y ) (p) (J, (p) dp dv ' 

Jo 1. J 0 


1 -2b A 

+ l BAv’) e 


W(p) <|> (p) 


X d »(l^) dp dV 


( 3 . 12 ) 


Eqnfi.(3*ii) and (3.12) simplify to 


fl -^a/\} 


-2bA 


^ 4 . => [ f yC^^o) +J^ ^ 


^ [ Vo ^C-v.) (e 


2^n/Vn -2b A„ -1 

0^ 0 . - O ’ 


, 4 . . 0 


)3 


(3.13) 
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and. 


If 

o 

o 

Vo 

" 2 

X(-V ) 

In [-. 0.:. ] 

X(v ) 

0^ 

Y(|J.) 



- .... . p 


2(1-0 

) X(-n) { 

f (c, v) = 

(1-c) 

g(c, v) x(-v) (v 

II 

> 

O 

[(1 - 

—1 ? 

c V tanh V ) + 


X(v) is the well-known X-function o.f Cass's theory and can 
be evaluated from the following eq.n. 

_ .'1 


X(-y) 


e3cp [~ ' g(c,v') (1 + ) In (v+v') dv] 

-^0 l - v '2 

( 3 * 16 ) 


Just like previous attempts on similar kinds of problems, 
oqna, (3.13) and (3.15) constitute the starting point of 
the present work. But instead of dealing the, two egtia-tions 
separately for numerical calculation, we combine them to .get 
a single non-homo geneous Bredholm equation for B_|_(v) 


B^(v) = f(o v)t“; y(uP + Y(|io) 

„ A -Sb/V 

+ I j ( V ) X(-v' ) e 

X V* (“'^^jr t h_^ (c,v, b)) dv*] (3.17) 


in which 
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li_j. ( c , V , b ) = 




V + \> 
0 


■-■2b7r“ 

0 


I'latheraatically, the slab albedo problem is defined exclu- 
sively by eon, (3.17). PJ’or a critical slab problem, 

1) in the absence of incident neutrons, = 0, and 

r 

T ( ) ~ ^ ^ ^ ^ M'q ) 

is identically zero, 

ii) because of the symmetry of the problem 


a_ . = 

0 + 0 - 

and A(v) ss A(-v) 

which implies b_ = 0, B__(v) =0. Since B_(v) is zero, 

B.(v) is replaced by B(v) to achieve simplicity in notation. 

T 

Hence eqn. (3*17) reduces to 

rl -2b/v' 

B(v) =s - ^ f(c,v)J X(-v')s 

X v' (-—i-,- + h (c, V , b) ) B(v’ ) dv' (3.18) 

V 4- 

where, 

f(o,v) = (o-l) g{o,v) X (-v) (v^ + k^) 

When 0 > 1, Vq ia imaginary and the fcernel of eqn. (3.18) 
can be further simplified by where is a real 

quantity. ' 
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Dofino, 


h ( o , V 5 b ) 


2z / V 

.‘I,° ° 

>-2b/v 

0 ° 

Vq*- V 


'2ir7v^ 

fj 0 0 

-2b/v 

+ e ^ 

2z^/ ik 
9. - - _ ° 

-2b/ik 
_ ^ 

it^-v 

ik + V 

0 


--2r-;iic7 — r5B}i¥ 

0 + c ' 


V, b) in rcGal and imaginary paris, wg got: 


Roal part of h_|_ (c, v? b) 


22+21 2z +2b 

{v (l + cos “■"“■g— — ) + sin . — — I. 


5 o 2z +2b 

+ kp {1 + cos (—9.^--.°) > 


b + z 


{ y + ]£ 

^ I ^ \ 0 

'*0 


( V ) O 


Imaginary part of b.^(c, v, b) = 0 
Thoroforo, 

1 


b + z 


h^(c, V , b) 


V + k^ 


{ V + k^ tan ■ —g* 


With this simplification, oq.n. (5.18) takes the form 

rl r^2b/v' 

° f {c,v) ; X(-v') e 

Jo 

V + tan -£ 


BCv) 



b+z. 



(5.19) 
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EqnB. (5.17) nnd (5.19) are the final form of Predholm eons 
we study in tho next section. 

5.2 Tho^ Inverse Eigenvalue Problem; 

Tho oqn.,(3.17) is rewritten in the general form of 
oigen value equation as, 


xE:j.(v) 

with, 


= p^(c,v, b) IC (v,v', b) Bjv') dv' (3.20) 

— Jo Z Z 

. -2b/v' 

b) =3 


X + li_|_ (c,v, b)) (3.21) 

and the inhomogeneous term i''_^ (c, v, b) is 

P_^(c,v, b) = f(c,v) Y(tAQ) ^ h^^(c,v,b)3 

(3.22) 


Eor critical problem (eqii. 3.19), the term E_|_(c,v,b) is absent 
and hence the eigenvalue problem is posed as a. homogeneous 
Eredholm equation, 


with. 


/'I 

XB(v) = \ K(v , v , b) B(v ' ) dV 
■.'o 


K{ V, v , b) « I 7(c,v) I(- V> ) e 


(3.23) 


.2b/v; 


b+z 


V + ICq tan ‘"jj* 


,2 + 

V + m 


V + V 


i-r ] 0 . 2 ^) 


X V ' [ 
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The integral operator genoratod by the kernel e:(v, v > , b) 


in oqn.(3.2'5) dofinoB a family of operators K(t) for difforent 
valuoo of the oigonvalue parameter b. Tlius the eigenvalue 


X (b) will bo dependent on b analytically i.e. as the ' 
parameter b changes, the operator !K(b) and its associated 
eigenvalue and eigenfunction change correspondingly. In the 


analysis of the ^.'igonvO.uo oejuation, it is determined how 
tho oigonvaluo X(b) and eigenfunction B(v) change with the 
operator !l<((b) i,o. v/ith changes in the parameter b. For 
a fixed value of tho paranetor b, the analysis reduces to 


that of usual. Fredholm equation and thus the solution of 
the oigonvaluo problem is straightforward. But, as we liave 
soon, tho problom arising in tho present study is the 
inverse of tho eigenvalue pi'oblem. In such a problem, the 
oigonvalu.: \ i.^ spc-cifiod and tho parameter b (i.e. tho 
operator ) which produces this eigenvalue has to be detor- 
miriod, Analysis of this Icind of problem, unlike the 
oigonvaluo problom, is indeed not simple and straightforward. 

5*3 Tho Critical Pro blom; 


In tho present analysis, the critical slab problem 
is defined by putting X (b) = 1 in equation (3.25) i.e. 

eigonvolue is equal to 3.. The parameter h corresponding to 
^ (b) a 1 is the oritical half-thickness and is the desired 
quantity of interest. The crux of the present analysis 
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coiifjis'ts in establishing the existence of such a b and 
SG'.'ks to determine this parameter. By putting X = i in 

GC^n, ( u ♦ 23 ) j txi^.' corresponding e^n, for the critical 
problem is obtained as, ■ 

rl 

B(’^) =\ K(v, v',b) B(v')dv' (3.25) 

**'0 X 

where K(v , v, h) is given by equation (3.24). 

To proceed with the analysis of eqn.(3.25), we need the 
following tliGorom due to Crohborg [4Q]. 

G ohb org Thq o r cm ; 

Lot iK('b) be an analytic operator va-lued function in 
an open connected set G, whose values are compact operators 
for b e G on the Banach space X. Then for any p 51^ 0, 

one of tho two possibilities mu.st hold (a) for every b s G, 

\i is an eigenvalue of 'K(b) or, (b) except for a discrete 
net of values b.,^ s G, the operator pH - IK (b) ,has a 
bounded every where defined inverse, while (u, il - !K(b)) 
has a pole at each of the points b^. t/e assume that the 
family of operators K(b) generated by the kernel K(v, v», b) 
acts, on an Hilbert space 1) of continuous positive 

functions. The spectrum of the operator il<( will, in general, 
consists of a continuous and residual part in addition to 
the isolated point spectrum which need not be finite. 

However, the compactness of the operator elimina^tes 
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the continuouB and residual spectra., and ensures point 
spcc'crura only. [T'lie into.;^ral operator iKCvjv', b) acting 


in 

l}^2 ( t , 

1) 

will be 

comp 0,1 

K( 

V , v , b) 

is a.n 'i 

-P kern 


1 

[•1 

1 .14(v ,, v' . 

b) ( 


0 J 

0 



Sine G V , 

V * 

e(0, 1), 

, it is 

be 

satisfied 

if tjic 

kernel 


(5.26) 


This is true in each of tho sogments. 


0 


and 


b + z 

0 , % 


<..^12. < (n-i-2) — Ih- n = 1,3,5,.,. 


In tliO derivation of ecin. (3.23), the slab thickness b has 

boon assumed to have positive real values. Mathematically 

also, negative values of b are not allowed. This is 

-2 b/v' 

because tho factor c in the kernel i:(v,v', b) in the 

neighbourhood of V =0, for negative b m.akes the kernel 

unbounded and thus oqn. (3.26) is no longer satisfied. 

Summarizing the results, we find that lK(b) is an ana-lytrc 

operator valued function in the open-connected’ ser G- either 

in the rang© 0 <_ • < ii:/2 or ^ < (n+2) 

0 o e 

n'=s 1,3, 5, ..whose values are compact operators for b e on 

trj2(0, 1), It can now bo concluded that the conditions for 
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applicability of 


Gohbere tboorem are satisfied 


in the present 


case. The 


b+i 

0 < -X 


0 


thoorc.ni is to bo applied in the range of 

% 

-2 ior our purpose. According to the theorem 


ono of the two possibilities must hold a) for every b e G, 

1 is the eigenvalue of the operator eqn. (5.25) B = E(b) B 
or b) except for a discrete set of bj^ s G, the operator 
f'— 5?](b) has a bounded evory“Wherc defined inverse^ while 
( f— {K(b)) has a polo at each of the points b|^, The 
operator 4K:(b) is obviously a bounded operator and hence 
from cqn. (3.25), 


It B II = IebII ^ K K )1 liB ti (3.27) 

For a nontrivial B to exist, it follows from eqn. ( 3 . 27 ) that 

ii eI' 2 1 

/I /I 

i.e. I > lK(v,v', b) r dv dv' > 1 (3.28) 

J ovo 

is the necessary condition for nontrivial solution of B. 
However, for b =0, it is seen from numerical calculation . 
that K(v,v', b) is alway.-) less than 1. Table 3.1 shows the 
values of K(v, v ', b) for ■ different values of c, and b = 0. 
Henco, for b = 0 1K|| is always less than 1 and it is 
concludsd,! cannot be an eigenvalue of eqn. (3.25) for b = 0. 
The first possibility (a) thus cannot be applicable and 
consoquently possibility (b) must hold. This proves the 
existence of atleast one bj^, (i.e. ihe critical thickness) 
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Tablo ! 


O.G 

-0,2689 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


).l(a); K(v, v' j b) as a function of 



V , V for 

c =1.1, 

b = 0. 


0.2 

0.4 

0.6 

0.8 

1.0 

-0.3697 

-0.2894 

-0.2387 

-0.2005 

-0.1674 

-0,1211 

-0.1241 

-0.1128 

-0.0988 

-0.0837 

-0.0562 

-0.0634 

-0.0600 

-0.0531 

-0.0447 

-0.0234 

-0.0333 

■ •0.0313 

-0.0279 

- 0.0227 

-0.0133 

-0.015S 

-0.0149 

-0.0127 

-0.0037 

-0.0000 

-0.0000 

-0.0000 

0.0000 

0.0000 
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l.'iolo ,),l(o); , V , 1) a,s a function of v, v ' 

for c = 1.5 and b = 0. 


V 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 5 

i 

0.0 

- 0.3150 

- 0.3830 

- 0.2580 

- 0.1797 

- 0.1195 

1 

- 0.0740 j 

CM 

• 

O 

0,0000 

- 0.1000 

- 0.0782 

- 0.0478 

- 0.0185 

1 

0.0068 1 

0 . 4 

0.0000 

- 0.0326 

- 0.0213 

C 

o 

o 

0.0147 

0.0313 ! 

0.6 

0,0000 

- 0.0095 

- 0.0017 

O.OIOG 

0.0218 

0.0329 . 1 

0.8 

0.0000 

- 0.0018 

0.0035 

0.0104 

0.0172 

0.0237 j 

1 

1.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 1 



I 

I 
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Tablo 

3 . 3. ( 0 ) : 

K(v, 

for 0 = 

, b) as a 

1.6, b = 

function 

0. 

of V ' 


\ I 

X 

V 

0.00 

0.2 

0.4 

0 . o 

0.8 

1.0 

0.0 

-0.3773 

-0.3650 

-0.1702 

-0.0436 

0.0504 

0.1233 

0.2 

0.0000 

-0.0483 

0.0147 

0.0775 

0.1297 

0.1744 

0.4 

0.0000 

0.0054 

0.0427 

0.0811 

0.1135 

0.1423 

0.6 

0.0000 

0.0122 

0.0349 

0.0582 

0.0779 

0.0958 

O 

CO 

0.0000 

0.0087 

0.0208 

0.0331 

0.0437 

0.0534 

1.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 






Table; K(v, V , b) as a function of v, v*, 

for 0 = 2 and b = 0. 


V \ 

0,0 

0.2 

0.4 

0.6 

0,8 

0.0 

- 0.4700 - 

- 0.3016 

0.0129 

0.2179 

0.3680 

0.2 

0,0000 

0.0402 

0.1668 

0.2746 

0.3623 

0.4 

0.0000 

0.0513 

0.1181 

0.1769 

0.2264 

0.6 

0.0000 

0.0317 

0.0664 

0.0976 

0.1245 

0.0 

0.0000 

0.0163 

0.0331 

0.0435 

0.0620 

1.0 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 


1.0 


0.4822 

0.4321 

0.2668 

0.1468 

0-0735 

0.0000 
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in the ranse of b e (0, | for vhloh 1 le the eigen- 

value of oc-n.(3.J‘5) end thus the criticality problem has a 
non— bri vial soluti- n. Tlio outcome of the results of Oobh erg 
thGorom loads to an important conclusion. It is to he 
noted that '2 ^q“^q tlie end-point critical half— thickness 
^ond-point present analysis it is proved that the 

exact critical half-thiclmess h. is atleast less than 
end-point • actually true as evident from exact 

critical thickness calculation by Sj^y-method. ) 

5 * 4 Kumoric al, Qa3 .ci^atioii 0 f Oritical Thickne s s by 

^octral Radius Methods 

As has boon mentioned in Oliapter I, the present 
approach of calculation of critical thickness is by no means 
a competitor to the fast converging Neumann series solution 
of Mitsie [ 15 ] » "but this method may, however, be addressed 
to such situo.tions whore numerical calculation of the 
inverse' para,metor is not otherwise possible or in cases 
where the Neumonn Series solution converges slowly, ¥e 
start with the concopt of spectral radius for a linear 
operator (K, It is proved in the theory of linear operators 
that tho limit 

r » lim 

n+® 

exists and is finite. The number r is called the spectral 
radius of the boundod linear operator 
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Th.0 sp-jd/Tfil 3?5idius T is ch.cirac'tGrized ly "blie facfc "biiai; "bh.© 
inecuality 

U I r 

implies that the operator ( iK - ^ !] )“^ exists and is 
hounded. In particular, if K is compact, then r = [x 1 , 

o ‘ ^ 

where Xq is an eigenvalue of E with maximum absolute 
value# In this analysis, the operator K as defined hy 
©qn# ( 3 * 25 ) ±3 a compact operator and hence the spectral 
radius r(h) is equal to the eigenvalue 1 i.e., 

r(h) = 1 . 

Itirther if tho operator K is a positive operator, fairly 
accurate ostimatos of houndo on spectral radius can he 
ohtainod. The oporator K is positive if it leaves a cone 
C invariant. A sot C E where S is a real Banach space 
is called a cone if the following conditions arc satisfied 
a) if u, V E C , then au + pv e t for all a, p ^ 0 h) if 
u e c and -u e C then u s: 0, The set of non-negative 
functions in a function space and the set of vectors with 
non— negative components are examples of cones. The criterion 
for positivity of an operator K then can he expressed 
hy saying that if x e fh then E x e <c . In case of an, 
integral oporator acting on the cone of non-negative conbi- 
nuous functions (as in our case) , the positivity of the 
operator is mot if the kernel is positive i.e, , 

K(v, v» ) > 0 
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for all 

V, V in their 

domain of definition n 

on the set 

of positive ron.1 numbers 

. Por such an integral 

operator. 

ii: can be. 

) nhovrn [ 49 ] 




r 

<. max 

ve n 

r 

■J 

fi 

r 

ii. ( V , V ) d V 


and , 



(3.29) 

r 

> min 

ve£2 

i 

n 

J.v(v , V * ) dv * 



I'h© coarse ostimatos of spectral radius can be further modi- 
fied as, • 

Jjz(v, v') Z( v' , t) dv* dr 

\ n 

^ tr t • N j . • 


r < max 
V en 


K(v , V ) dv 


and * 


'I j Z(v, v') Z(v', t) dv' dT 


r > min iLi? — — 

j Z(v, v') dv' 
£2 


(5.30) 


Prom eqns. (3.29) and (3.30) the actual spectral radius can 
be approximated as the average 


r * 


1 [ max 

2 V e 


and 


r * 


w [ max 
veQ 


+ rain 


( I\( V, V ' ) d v' -I- min 5 K(v,v ) dv ] 

vef2 n 




jJzCv, v') K(v', t) dV dT 


(3.31a) 


J 

£2 n 


£2 

j^K(v, v')K(v','t) dv' dT 


iK(v, v')dv' 


] (3.31b) 
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Til© 0 pud 6 si>pppo3ciiii8."tioii of sp©©^!?!! rsidius givsu 
by oqn. (5.51a) s however, gives a fairly accurate estimate 
of the la'-gent eigenvalue of a: Fredholm operator as the 
followin’^;; examples show. The eigenvalue obtained by 
eqn, (3.51a) io denoted as "x - 


a) 


'^(v, v’) 


o'*-'-" J 1,, v’ e (0, 1) 


h) 


^ (e^-l) = X (exact value) 

;C(Vj V ) = V + V ' } V , V ' £( 0 , 1 ) 


X e 

I 


c) 


X s 1.00, X = 1.077 
K(v, V ) =s ^ v'(2«-v)} V, v' e (O, 1) 


0.575 X = 0.333 


The estimates of spectr"'! radius as given by eqn. (3.31a) 
will now be applied for the criticality eqn. (3.25) to 
obtain an approximate numerical estimate of critical 
thichness. 


The eqn. (3.25) is i-ewritten here for the purpose of 
analysis, , 

B(v) mi K(v, v', b) B(v') dv ^^^^^^^^^^^^^^^^(3-25) 

Jo 


With, 


K( V, V #b) a 2 

X V C 


° f (c,v ) X (-v' ) e 

b+z 

v+ k tan -p- 
o 


.2b/v' 


0 


2 + 

V + k^ 


v-fv 


(3.24) 
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Unfortunately, the kerne], given by eqn. (3.24) is not ' 
positive for all v' s (o, i). positivity of , the 
kernel will be met if 


V + k„ tan 


b+55 

( 


2 1 2 
V + k 

0 


v+v 


i.o. if 


..2 I 1.2 

V + ilu_ 

0 

“"Ij+Z 
V + k^ tan -- 1 ^— 
0 


b+2 


k‘ 


V k^ tan 


b-i-2 


V + k^ tan 


Since the right hand side has a maximum value at v = 0 for 
V c (0, 1) the positivity criterion reduces to 

k 


_o_ 

b+2. 


tan (-ip— ) 


Thus the kernel K( v, v'^ b) (and hence the operator E) is 
porjitivo for all b in the range of v e(0, l) and v’ belonging 
to the sub~set (vj^(b), 1) of (0,1) where 

■ k . ■ 

vjih) = 

b+2 

tan 


(3.52) 
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Howovor, 

) 

of V, V 


for 3'redb.olni operator^ the domain of definition 
Bhould be the same. Hence both v, v ' are allowed 


only in. cb.<i iiit0.i.'V.il (vy^(b)j l)j which is obviously a 


conservative estimate. (In appendix, it has been shown that 
further lowering of Vj^(b) is, indeed, possible.) Prom 
eqn. (3.3^?), it follows, 



By putting raaximum posoible value of = 1, eqn. (3.33) 
gives a lower bound for critical half-thickness as 


\ tan“^(k^) - (3.34) 

wliaxfias tha upper bound has already been derived from 
Gohberg theorem as, 

”0 = \ 


At the end-point critical half-thickness, it is observed 
from eqn. (3.32) that = 0. Since the end- 

point result is fairly accurate, it can be reasonably assumed 
that \^,(b) will not be substantially different from zero 
when b is the exact critics, 1 half-thicloiess. This fact 
wi.ll be reflected by numerical calculation to follow in the 
next section. 

In order to ensure positivity of the operator K, 

we had to slightly deviate from the original problem, x.e,. 


i 
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¥0 do not nolvo the exact prohlem aa posed hy eqn. (3. 25) Put 
tlio following approxinato ono 


where . 


B(v) 






1 




K(v, v', b) B(v’) dv' 


Prom oqn* (3.31a), it follows 

rCh) = 




Whore , 


and 


r (b ) 

raax^ 




raax '• K(v, v , b) dv 

ve (0,1) 1 '^j^(b) 


fl 

min ( IC(v , V , b) dv 

V e(0 ,l) Jvj (b) 


0.'he critical hnlf-thicknesa is given by the eqn. 


>(b) 


The inte^-ral 


1 


'vj^^(b) 


(3.35) 


(3.36) 


Xiv, v', b) dv' can be evaluated as 


Kv, b) 


fl I I 

\ .'f ( V , V , b ) dv 


-'v^(b) 


b+z 


I f (c, v) 


y -I- k tan --g" 

^ r 


2 , 2 
V + h 


ipb) 


E ‘f (c,v ) 


(3-37) 
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wlioro, 

r 1 - 2b / V 

I-{(b) .. v’ e dv’ 

(""1 ' — 2b / V ' 

and, Io(v, b) = .„^X(_v')e dv ’ 

-‘^.j.(b) v+v 

The integrals l 2 _Cb) and l 2 (^» b) can be obtained numerically 
with X(~v ) given by eon, (3.16). However, to get arousid 
the ■time-consuming numerical integration, I^Cb) and 
IpCv, b) are approximated in the present analysis as 

-2b 1- v^(b) 

l 3 _(b) « X(-l) e ( ) 

IpCv, b) «. X(-l) [1 - Nv(b) - V In ] 

V + V j^(b) 


Thus, I(v, b) is approximated as 


b+z 


V + k tan 


I(v, b) 


0 

7 


.2b 


IT*' 


« e - ~ Z ( - 1 ) f ( c , V ) { — -72 


2(v^ + k^) 


X (.1- M?{13)) - U 

v + v ^ vd ; 

d.i • 

The constants k^, g(c, v ) needed for evaluation of I (■w, b) 

in eqn, (3,38) are taken from Case, de Hoffmann and Plac zee k [23] 
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Eqn. (3.35) is infact a transcendental equation 
in b B.n^ the bir^ection method is employed for evaluation 
of b. A vraue of b = b_^ in first chosen in the range of 
bj^< b < bg and I(v), b) is calculated from eqn. (3.38) 
for 21 points of v in tne interval (^jj> !)• The maximum 
of these values give while the minimum corresponds to 

%in‘ Bpectral radius is then calculated from eqn. (3.35). 

If r(b-j_) is equal to 1 within a specified error limit (in the 
present calculation this lj.mit is taken as e = | r(b)-lt 
= lO”^)* then this choice of b = b^ is the desired critical 
hcilf~thicknees b^. If tiK; convergence criterion is not met 
i.e. ^(b^^) It then another value of b = b 2 is chosen 
such that r(b 2 ) - 1 is of opposite sign to r(b^) - 1. The 
bisection method is then continued to achieve the desired 
convergence in spoctral is.dius. The corresponding value of 
b is the critical half-thickness b^. Numerical results of 
this calculation alongt'^itli and values are shown in 
Tabic 3*?.. values at b = b^ for different c are also 

given in the table. It in observed from numerical values 
of that the approximate problem does not differ significantly 

from the original one for the range of practical values of 
c ( jj 1,6), Considering the crudity of bounds of spectral 
radius used here, the results obtained by the present meuhod 
can be considered to be fairly accurate. The error in critical 
thickness increases with increasing value of c. This fact is 
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Table 

3.2; Oiltical half-' 
moan froo path 

thickness b 

• 

Q in units 

of 







0 

PrsBsnt 

error) 

i'fo error) 

P j 

(/> error) 

Exact 


1.02 

5.6655 

(0,000) 

5.7519 

(-1.525) 

5.6711 

(-0.099) 

5.6655 

0.0000 

1.09 

3,3002 

(0.000) 

3.4034 

(-3.438) 

3.3066 

(-0.193) 

3.3002 

0.0000 

1.1 

2,1125 

(0.043) 

2.2287 

(-5.454) 

2.1213 

(-0.376) 

2.1134 

0.0008 

1.2 

1,285^^ 

(0.295) 

1.4125 

(-9.555) 

1.3020 

(-0.985) 

1.2893 

0.0043 

1.4 

0,7231 

(1.833) 

0.8581 

(-15.54) 

0.7577 

(-2,860) 

0.7366 

0.0137 

1,6 

0.4927 

(3.769) 

0.6269 

(-22.44) 

0.3394 

(-5.151) 

0.5120 

0.0225 

1.8 

0.3634 

(6.509) 

0.4965 

(-27.73) 

0.4182 

(-7.592) 

0.3887 

0.0298 

2,0 

0,2805 

(9.749) 

0.4121 

(-32.59) 

0.3420 

(-10.05) 

0.3108 

0.0359 
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not 3uri: rising. With incroasrng value of c, deviates 
more from '/.oro and naturally the approximate problem is 
not that c.lo.’o to tho original one as compared to lower 
valaeo of c. 'lowevor, the critical thickness significantly 
iraprovea if intogrationo in eqn. (3»37) are perfomed 
numerically instead of using eqn. (3.38). In this case, 
the error for c = 2.0 reducos to 7.876 percent from 9.749 per 
cent. Nuraorical integration is, however, not recommended 
beoauso of much greater computational time. 

3.5 yh.9. .PAfiA, iy -.!bp.<3„o„ .Prphlpm 

For tho analyeis of the slab albedo problem we 
consider the oqn, (3.20) with X = 1 

fl 

(y) « (c,v, b) -j E:^(v, v', b) B^(v') dv' (3.39) 

The kornol K,.(v, v’, b) is givan by eqn. (3.21). PorB_^(v), 
K^(v, v', b) is always nositivo for v, v' e (0, 1) and 
for all b. Thus the homogeneous eqn. 

A 

B^{v) « ~j Kj.(v, v', b) B.|.(v’) dv' 

oamiot have a nontrivial solution in the cone of positive 
functions as the eigenvalue in this case is -1. Since B_j_('') 
is a linear combination of A(v) and A(-v) as seen from 
oqn, (3.10), A(v) and A(-v) are also non-existent and thus 
B„{v) of the homogeneous equation also does not exist. By 
Fredholm alternative theorem [50], it can then be concluded 
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tlmt 01 inhomogonooun oqn. (3.39) have unique non- 


trivial :;n:i uM(Viij:;i. 


.rjOo';:ing this from another angle, this 


can also r.hown o.a follows, 3__(v) satisfies the following 


equation, 

('1 

-M K (v, v', h) B (v') dv' (3.4Q) 
•0 

whoro the kornol K_^(v, v', h) of eqn. (3.21) is redefined as 

, I V o -2b/v' 

K„(v, V , h) =1 f(c,v) I(™v') e 


X v' H- h (c,v,h)) 

v+v 

¥c consider tho governing aqn, for B_(v ) only as we find 
IC_^( V, v' , "b) > K^(v, V *, To) and any inference drawn on 
B_^(v) in this analysis trill hold for B_j_(v) also. To show 
that B__(v) has a unique non-trivial solution, we proceed 
as follows: 

An equation of the form 
XX = IK X + f 
has a unique solution 

X a (X I - ^ 

if I ^ I > spectral radius 

From oqn, (3.29), the upper hound for spectral radius r can 
be estimated in our case as, 


(3.41) 

(3.42) 

(5.43) 
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= max I^(v, b) 

ve(0,l) ^ 

The right hand side I^( v, b) is calculated numerically for 
different values of c <1 and for different values of b 
for 21 points of v in (0, 1) starting from v = 0, at an 
interval of ,05. The results are shown in Table 5.5. 
Maximum value is always found to be less than 1 and hence 
the spectral radius is less than 1. Prom eqn, (5.40), it 
follows Ia| = 1 and so, condition (5.45) is satisfied. 
Hence, eqn. (5.40) has a unique solution for B_(v ). In 
other words, the slab albedo problem has a unique solution 
for all c < 1 and for all b. The restriction relating b 





Table 3.5(a) 


Table 3 

.3(a) I^(v, 

b = 0 

b) as a function .of 
. 5 and c = 0.9? 0.6, 

V for 

0.2. 



c 


V *“* 

0.9 

0.6 

0.2 

0..00 

0,0848 

0.0649 

0.0256 

0.05 

0,0676 

0.0536 

0.0219 

0.10 

0,0578 

0.047? 

0.0197 

0.15 

0,0507 

0,0425 

0.0181 

0.20 

0,045?- 

0.0389 

0.0168 

0.25 

0,0407 

0,0360 

0.0153 

O 

• 

o 

0.0369 

0.0337 

0.0149 

0.35 

0.0337 

0.0317 

0.014? 

0.40 

0.0309 

0.0301 

0.0136 

0.45 

0.0284 

0.0288 

0.0131 

0,50 

0.0265 

0.0276 

0,01^27 

0,55 

0,0241 

0,0267 

0,0124 

0.60 

0,0222 

0,0259 

0.0121 

0,65 

0,0 204 

0.0253 

0.0120 

0.70 

0,0137 

0.0248 

0.0119 

0.75 

0.0169 

0,0244 

0.0120 

0.80 

0,0151 

0,0241 

0,0122 

0.85 

0,013? 

0.0238 

0,0127 

0.90 

0.0110 

0.0232 

0.0136 

0.95 

0.0082 

0.0212 

0,0159 

1.00 

0.0000 

0.0000 

0,0000 
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Table 

3.3(b): I^(v, 

and c 

b) as a function 
= 0.9, 0.6, 0.2. 

of V for b = 1 

V 


e 



0.9 

0.6 

0.2 

8 

o 

0.0293 

0.0234 

0.0094 

0;05 

0.0232 

0.0193 

0.0081 

0.10 

0.0197 

0.0170 

0.0073 

0.15 

0,0172 

0.0153 

0,0067 

O 

C\J 

• 

O 

0.0153 

0.0140 

0,0062 

0..25 

0,0137 

0,0130 

0.0058 

o 

* 

o 

0,. 0124 

0.0121 

0.0055 

0.35 

0,0113 

0.0114 

0.0052 

0.40 

0,0103 

0.0108 

0.0050 

0.45 

0,0095 

0.0103 

0,0048 

0.50 

0.0088 

0.0099 

0.0047 

0.55 

0.0080 

0.0096 

0.0046 

0.60 

0.0074 

0.0093 

0.0045 

0.65 

0.0068 

0,0091 

0.0044 

0.70 

0.0062 

0.0089 

0.0044 

0.75 

0.0056 

0,0087 

0,0044 

0.80 

0,0050 

0,0086 

0.0045 

0,85 

0.0043 

0.0085 

0.0047 

0,90 

0,0036 

0.0083 

0.0050 

0.95 

0.0027 

0.0076 

0.005s 

1.00 

0.0000 

0,0000 

0.0000 
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and c as obtained by McCormick and Mendelson [17] is 
only necessar 3 ' for the Senmann Series solution to converge 
and not for a solution to exist. 


CHA.PTER IV 


MODIPICA.TIOI AND STJGG-ISTIONS 


Intr oduct ion; 

In this chapter, an attempt has been made to rectify 
some of the shortcomings as appeared in the development of 
Chapter II, It is concluded in Chapter II that the 
formulalism based on source-free Milne problem does not 
necessarily improve the leakage for the constant source 
Milne problem for strong absorption cases. The reason for 
tliis is attributed to the fact that the extrapolation 
distance for constant source problem is widely different 
from that of source-free Milne problem for low values of 
c. This necessitates a separate formulation of 
approximation based on constant source Milne problem in 
those cases where the source-free Milne problem cannot 
bo considered as the standard reference problem. Here, 
an expression for the extrapolation distance for constant 
source problem has been derived by application of Case's 
theory, TP^ formulation based on constant source problem 
is then constiu.cted by introducing a parameter in the 
Marshak free-surface condition in a similar manner as in 
Chapter II. The parameter is determined by equating 
the leakage from the half -space in TP^^ approximation to the 
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exact leakage, Tlie formalism car be easily extended for 
high-or order approximations. However, it has not been 
undertaken in the present work. The purpose of the present 
work is simply to demonstrate that the source-free Milne 
problem should not be taken as the only standard reference 
problem for construction of the TP^j equations. The scheme 
can be applied for any practical problem which can be appro- 
ximately represented by a constant source type problem such 
as computation of thermal worth of a control element in an 
assembly. In this chapter, the escape probability for a 
finite slab containing a source has been calculated by the 
present method and compared with the results of conventional 
Pi approximation and exact analysis. It has been demonstrated 
that the present formula,tion improves the results whereas 
the TPj^ approximation based on source-free Milne problem 
fails for strong absorption cases. 


4.1 Anal ysi s ; 


For the constant source half -space problem, the 


extrapolation distance is obtained by equating the asjnapto- 


tic flux to zero as follows [5] 


2qo 

TTPcT 


+ a 


04 " 




where. 


(4.1) 


■% 
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cl 




0 + 


(1-c) c 


“n 2(vJ 


Is a uniform isotropic source 
Vq is given by the positive root of 


CV V+l 

jLiI m 1 

d y 


We already have. 


X(v^) = 


z /V r 
' \/ 


v; (1^.0) 


2v2 (i„c) (v2 _ 1) 


(4.2) 


(4.3) 



is then ohtaine(3 from eqns* 


cz 


cs 

"0 


c V 

In { 


(4.1), (4.2) 
[^0 

2(l-c) (v^ ~ 


and 

+ 13 
1 ) ' 


(4.3) as 

> (4.4) 


Table '4.1 compares the values of as obtained from 

oq.n, (4.4) with z^. It is found that for low values of c, 
c s 

z^^ is widely different from z^. Thus , TP^-approximation 
normalized to exact asjj'mptotic flux of sourco-free Milne 
problem (i.e, exact z^) may lead to significant errors for 
constant source type problem in case of strong absorption. 
This calls for a separate formulation for this class of 
problems. However, unlike TP^ theory developed in Chapter II, 
tho present formalism is not normalized to exact aS3nnptotic 


Ta'blc 4.1; Constant-Source Milne problem 

c s 

extrapolation distance and 
source-free Milne problem 
extrapolation distance z^. 


e 

cs 

cz 

0 

CZ 

0 

0.0 

2.000 

1.000 

0.1 

1.4097 

0.8559 

0.2 

1.1162 

0.7851 

0.3 

0.9439 

0.7491 

0.4 

0.8418 ' 

0.7305 

0.5 

0.7792 

0.7205 

0,6 

0.7401 

0.7155 

0.7 

0.7158 

0.7127 

0.3 

0.7023 

0.7113 

•0.9 

0.6977 

0.7106 

1.0 

0.7104 

0.7104 
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fluz. The- reason for this is the following. The use of 
exact asymptotic flux in constant source problem yields an 
infinite linear extrapolation distance for the purely 
absorption case [51]. It has been found [51] that the use 
of this linear extiapolation distance when used in asymptotic 
diffusion theory would yield zero control rod worth. In 
order to circumvent this unrealistic situation, the present 
formulation is normalized to exact leakage from an half-space 
and not to the exact asymxototic flux. It is to be noted 
that this problem does not arise for . source-free Milne 
problem in which case the extrapolation distance has a 
finite value as c -^0. 


TP^^ flux for the constant source Milne problem is 
given as, 


’f' (x, p) 


(|) fn) 



(4.5) 


The leakage from the half-space is obtained from the above 
eon. as , 

'1 




p 4-(0, p) dp = a^^^ g^(v^) (4.6) 


The Ifershak free-surface condition is modified by introducing 
a parameter following manner. 


1 (p+E^) ’l^(0> p) <3p 

J 0 


(4.7) 


0 
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By- use of eqn. (4.5), the expression for is 
obtained from condition (4.7) as 


^1 


- V 
2 ( 1 - 0 ) 


P™ { Vr 

■ “^0+ ” ^ ^O “IT 

o 


(1) ^''o 

TT-^T " ^0+ “T- 
( 1 ) . 


(4.8) 


"f a 


In 


'■o-l 


The coefficient is obtained by forcing the TP^^ leakage 

equal to exact leakage as 


(Tl) _ ^ 

exact 


2q 

- 1(c)] 

C 0 V / 


(4.9) 


By use of eqn. (4.5), a^^'^ is obtained from eqn. (4.9) as 

2C! 


0 + 


- !(=)] 


( 4 . 10 ) 


Substituting a^^^ from eqn. (4.10) into eqn. (4.8) e. 


simplifies to 


( V„-I(C)1 {V^ _1,. 1 


e-,(c) = 


(4.11) 


{v^ - I(c)> In - 1 

o 


The function 1(c) is tabulated by Case, de Hoffmann and 
Placzeck [23]. The limiting values of e^(c) for pure^^ 
scattering case is obtained by a limiting procedure to be 


eqd) 


2_- = 

3(21(1) -1) 


- 0.1039 
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In. ondor "bo avoid a loiigiliy liinii;ing procedure for pure 
absorption case, £ 2 _( 0 ) is simply evaluated by extrapolation. 
I'liis is found to be equal to —0,0035. Due to extreme 
insensitiveness of value with c for bigh absorptions, 
it is recommended that £^^( 0 ) should be evaluated by simple 
linear interpolation for values of c between 0 and 0,2, 
Uumerical values of for different values of c are given 
in Table 4.2. For any problem involving freo-surface and 
representing a constant source-type problem , e^(c) as 
evaluated from eqn, (4.11) is used. The coefficient ai^^ 
can then be obtained from Marshak free-surface condition (4.7), 
The formalism is now applied for computing escape probabi- 
lity P for a finite slab. 

^*2 Escape Probability for a Pinite Slab; 


¥e consider a homogeneous slab of half-thickness 
b with a spatially constant source and the rest of the space 
considered as a purely absorbing medium. The angular density 
'P (x, p) for TP]_-approximation, with the origin of coordinate 
at the centre of the slab, is given by, 


4 (x, p) = 


rr-^ + < 


cv 


2(v^-p) 


■x/v. 


cv 2;/v^ 


2(v '!-p) 


(4.12) 
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Talle 4.2; £^(0) a,.s a function of c. 


c 


0.0 

-0.0035 

0.1 

-0.0074 

0.2 

-0.0113 

0.3 

-0.0164 

0.4 

-0.0228 

0.5 

-0.0292 

0.6 

-0.0368 

0.7 

-0.0445 

0.8 

-0.0544 

0.9 

-0.0679 

1.0 

-0.1039 
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IhG coGfficient is oltainod from the condition (4.7) as 

I (P- + s-j) IJ.) dfi = 0 

-^o 

By use of eqn. (4.12) for (-b, p.), we got 

i,(l) _ 


0 + 


^ 

b/'' V “b/v ,, J.T 

CV (I_c) { (v +£ ) o ° In --°y - (v -e )e ° In 
o o ± V — -L o ± • , 

0 '’o 

- 2 sinh(b/v^> (4.13) 

In oqn. (4.13), is used from eqn.(4.1l). The escape 
probability P is defined by, 


f 

(1-c) j <!> (x) dx 


p ^ 1 „ 


4Dq^ 


(4.14) 


The scalar flux b (x) from eqn. (4.12) is 


f 1 ^ 

(|)(x) = °~ + cosh(x/vJ 


(1-c) 


(4.15) 


By use of eqns. (4.14) and (4.15), P is given by 

^ 0 +^ sinh (^ ) g-, (v^) 


,(T1) 


o ^ ^1^ o 




(4.16) 


Substituting a^^^ from eqn. (4.13), oqn. (4.16) becomes 


(Tl) 


(1 + 2e^) Sinh (-J^) 
Vo 


-b/v^ 


+1 


cb {(Vo+Si) 6 ° In — JY -( V^l^ ® ‘ 


V 


- 2 sinh(b v^) ^ (4.17) 
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Convontional P^-approxLmation result is given ty, 




b { 


cosn 


';rdT> 


+ 2 


^1 




sinh (- 


b 

Tl7 


)> 


Table 4.3 compares the TP^^ " values of escape probability vxith 
Pj_v and exact values. TP^ values corresponding to source- 
free Milne problem are also given for comparison. It is 
found that the present formulation improves the result 
significantly over conventional P^-approzimation for all 
values of b and c where as the TPj^^-formulation based on 
source-free Milne problem fails absolutely for low values 
of c. The exact results for c = 0,6, 0.4 and 0,2 are taken 
from curves [52] and hence, they could be given only upto 
two decimal places. The other exact results are taken from 
rcforence [26], 

'^•3 S uggestions ; 


It has been observed that TP formulation based on 
source-free Milne problem fails in those cases where the 
transient is overwhelmingly predominant., This is because of 
the fact that the transient is still in the framework of 
I’u approximation and the improvement in the assnnptotic portion 
of the solution is offset by a corresponding worsening of 
tho transient part. The conventional Pjj approximation, 
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Table 

4.32 Escape probability from a slab, P, 
function of b and c. 

a-S a 

c b 

b ■ 

{fo error) 

based” on ” ~ 
constant source 
Milne problem 

error) 

based" 'on””' 
source-free 
Milne problem 

Exact 

•0.5 

0.5024 

(4.6l) 

0.8402 

(-2.59) 

0.8267 

(-4.16) 

0.8626 

0.9 2.0 

0. 5866 

0.5425 

» 0 . 5300 

0.5449 


(7.65) 

(-0.44) 

(-2.73) 


10.0 

0.1400 

(7.36) 

0.1305 

(0.07) 

0.1272 

(-2.45) 

0.1304 

0*5 

0.8223 

(8.34) 

0.7237 

(-4.65) 

0.6839 

(-9.89) 

0.7590 

0.8 2.0 

0.4233 

(10.72) 

0.3803 

(- 0 . 52 ) 

0.3547 

(-7.22) 

0.5823 

10.0 

0.0903 

(9.85) 

0.0823 

(0.12) 

0.0766 

(-6.81) 

0.0822 

0.5 

0.70 

0..56 

0.44 

0.61 

0.6 2.0 

0.28 

0.24 

0.19 

0.26 

10.0 

0.06 

0.05 

0.04 

0.05 

0.5 

0.6499 

0.5016 

0.3316 

0.5591 


(16.24) 

(-10.28) 

(-40.69) 


0.5 2.0 

0.2408 

(14.83) 

0.2082 

(-0,71) 

0.1346 

(-35.81) 

0 . 2097 

10.0 

, 0.0486 
(14.08) 

0.0427 

(10.23) 

0.0276 

(-35.21) 

0.0426 

0.5 

O.Sl 

0.45 

0.23 

0.52 

0.4 2.0 

0.21 

0.18 

0,09 

0.19 

10.0 

0.04 

0.04 

0.02 

0.04 

.. 0.5 

0.54 

0.37 

0.09 

0.45 

0.2 2.0 

0..17 

0.15 

0.04 

0.16 

10.0 

0.05 

0.03 

0.01 

0.05 


/ 
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being a total theory, may however predict the overall 
results fairly accuratoly although the contributions due 
to asymptotic and transient parts separately arc- far from 
the exact results, The improvement of the fpj^-foimulation 
for strong absorption cases, however, may not be 
possible in the framework of -approximation as’ has been 
discussed in Chapter II. A better transient flux repre- 
sentation containing few teims is thus called for. 

>Spectral radius method of inverse parameter estima- 
tion can be applied for some general problems in those 
cases where the inverse parameter is not otherwise easily 
available. 

Ho such attempt lir^s, however, been imdertaken in the 


present dissertation. 
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APPENDIX A 


CONSIRUGIION OP ) PUNOTIONS POR ANISOIROPIC 

SGAITERING 


In case of anisotfopic scattering tlie one-speed 
tianspoft eqtiation is, 

1 

I 

1=0 

a 


P 


JiX^.Pi + ^ ^ c (21+1) f-, P, (^) 

& X ^ -t + J- 


X f f (x, tJ,' ) Pj_ (p' ) <3p' 

J -1 


(Al) 


Where the scattering function is approximated by a finite 
sum of Legendre polynomials and L is the order of anisotropy. 
In a similar procedure as in isotropic case, we apply the 

ansatZj 


“XA 

ii« (x, p) = e <!) y(p) 


(A2) 


where. 


Ui\^) = I S^iv) \ip-) 


i^3) 


n=o 


and, obviously, 


iv 


) = f <!» v(p) 5*11 (p) 

J-1 


CA4) 




By use of Eqn. (A2) in (Al) 4'y(|j.) is obtained as, , 


T. 




C V 


2 (v “[I ) l(=o 


I (21+1) f3_ g^Cv) 


i (- 1 , 1 ) 


(A5) 


and , 




+ X ( V ) 6 (v — [.i ), V e(-l,l) 


(A6) 


Use of eqns. (A3) and (A6) gives g^(v ), when v e(-l, l), as 


I. (-1 Pi(^) ^ 

g^(v ) = J (21+1) f^ gT_(v) P j _ ^ 


1=0 


-1 (V-^l) 


+ X(v)p^(v) 


, V £ (-1, 1) 


By utilizing, [53] 

Pi(i-i ) P^(^t) 

^ -1 (^> -la) 

¥e get, when ^ e (-1, l) 


dp, = 




g. 


‘n 


(v) = CV P Q^(v) I (21+1) f,_ g^(v) P^(v) 


1=0 


+ X (v ) 


n > L 


(AS) 


and, 



Sn(v) = ^ {PQ^C'') J (2-c\^)f21+l) g^(v) Pj^tv) 

L 

+ (21+1) g^(v) PQ^(v)l 

+ X (v) Pj^(v), n < I (A9) 

The factor X(v) is defined from Bqn. (A9) with = 1, i. 


X(v) = 


1 


I. 

c vP ( 21+1 ) 

1=0 


fl g]_(v) Q^(v) 


when ^ (-1» 1) j x(v) = 0 and the asymptotic roots are 
obtained from the following transcendental eqn. 


1 

I 

1=0 


c I (a+1) f-. 


Qi(-) 


= 1 


(AlO) 


B. Modification of 

It has been shox-m in Chapter III, that the positivity 
of K( V , V ’ ^ To) is met if 


? b+z 

ko - v^o tan 

lo'+z" 

v+ k^ tan 


(Bl) 


This condition is satisfied for all values- of 
ic 

yi e (™ e(0, l) as has been shown in 


tan 


"'TT 



Chapter III, However, for the puip)ose of 
lower limit of has also "been changed to 

is obviously a conservative estimate. It 
here that if a new lower limit of ^ and 


analysis, the 
k 


h+z 

tan -‘i. 


r-- , This 


will be shown 
v* is chosen as 


L - 


2 


k 


b+z. 


2 tan 


the positivity criterion (Bl) is still satisfied, 
hand side of (31) has. then a maximum value at 
which is eq.ual to 


The right 




- 

o 


k 


b+z 


o 1 , 0 

-T+z - • “IE"' 

2 tan ° 

jX 


k 


b+z 

2 tan - 

x\. 

o 


b+z 

+ k^ tan 


b+z^ 

k tan --r- 
° '"o 


’b+z 

(1+ 2 tan^ - ' ■"-) 



b+z 

( 2 tan ' ) 

o o 



Thus, it is found that the condition (Bl) is met for all 
•values of v’ e expected that this new 

lower limit will induce a correction in 

critical thickness r esults* further lowering of may 
also he possihle. 

However, as no estimation of spectral radius for 
non-positive operators may he possihle hy existing 
theories, the positivity criterion has to he satisfied in 
any cas e , 




